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I review the history of HBT interferometry, since its discovery in the mid 50's, up to the recent 
developments and results from BNL/RHIC experiments. I focus the discussion on the contributions 
to the subject given by members of our Brazilian group. 
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I. INTRODUCTION 

I will discuss here the fascinating method invented 
decades ago, vifhich turned into a very active field of in- 
vestigation up to the present. This year, we are celebrat- 
ing the 50"* anniversary of the first publication of the 
phenomenon observed through this method. In this sec- 
tion, I will briefly tell the story about the phenomenon 
in radio-astronomy, the subsequent observation of a sim- 
ilar one outside its original realm, and many a posteriori 
developments in the field, up to the present. 



I.l HBT 



Hanbury-Brown, as a means to measuring stellar radii 
through the angle subtended by nearby stars, as seen 
from the Earth's surface. Before actually performing the 





FIG. 1: Aerial photo and illustration of the original HBT 
apparatus. They have been extracted from Ref/l-]. 



HBT interferometry, also known as two-identical- 
particle correlation, was idealized in the 1950s by Robert 




FIG. 2: Picture of the two telescopes used in the HBT exper- 
iments. The figure was extracted from Ref.l^. 

experiment, Hanbury-Brown invited Richard Q. Twiss 
to develop the mathematical theory of intensity inter- 
ference (second-order interference) 2]. A very interest- 
ing aspect of this experiment is that it was conceived 
by both physicists, who also built the apparatus them- 
selves, made the experiment in Narrabri, Australia, and 
finally, analyzed the data. Nowadays, the experiments 
doing HBT at the RHIC/BNL accelerator have hundreds 
of participants. We could briefly summarize the experi- 
ment by informing that it consisted of two mirrors, each 
one focusing the light from a star onto a photo-multiplier 
tube. An essential ingredient of the device was the cor- 
relator, i.e., an electronic circuit that received the signals 
from both mirrors and multiplied them. As Hanbury- 
Brown himself described it, they " ... collected light 
as rain in a bucket ... ", there was no need to form 
a conventional image: the (paraboloidal) telescopes used 
for radio-astronomy would be enough, but with light- 
reflecting surfaces. The necessary precision of the sur- 
faces was governed by maximum permissible field of view. 
The draw-back they had to face in the first years was the 
skepticism of the community about the correctness of the 
results. Some scientists considered that the observation 
could not be real because it would violate Quantum Me- 
chanics. In reality, in 1956, helped by Purcell 3], they 
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managed to show that it was the other way round: not 
only the phenomenon existed, but it also followed from 
the fact that photons tended to arrive in pairs at the two 
correlators, as a consequence of Bose-Einstein statistics. 
A very interesting review about these early years was 
written by Gerson GoldhaberQ, one of the experimen- 
talists responsible for discovering the identical particle 
correlation in the opposite realm of HBT: the niicrocos- 
mos of high energy collisions. 

1.2 GGLP 

In 1959, Goldhaber, Goldhaber, Lee and Pais per- 
formed an experiment at the Bevalac/LBL, in Berkeley, 
CA, USA, aiming at the discovery of the p° resonance 4]. 
In the experiment, they considered pp collisions, at 1.05 
GeV/c. They were searching for the resonance by means 
of the decay p° tt+tt^, by measuring the unlike pair, 
TT'^TT~ , mass-distribution and comparing it with the ones 
for like pairs, tt^tt*. Afterwards, they concluded that 
there was not enough statistics for establishing the exis- 
tence of p". Nevertheless, they observed an unexpected 
angular correlation among identical pions! Later, in 1960, 
they successfully reproduced the empirical angular distri- 
bution by a detailed multi-7r phase-space calculation us- 
ing symmetrized wave functions for LIKE particles. Be- 
ing so, they concluded the effect was a consequence of 
the Bose-Einstein nature of tt'^tt'^ and 7r~7r~. They were 
not aware of the experiment Hanbury-Brown and Twiss 
had performed previously. Thus, they had discovered, by 
chance, the counterpart of the HBT effect in high energy 
collisions. They parameterized the observed correlation 
as: 

^1 ^2 2 / 2 2\ 2 

C(g2) = l + e-'S'- =i + e(9o-q)'^ (1) 
- -q^ = -{ki - k2f = Ml^ - (mi + m-if . 

The Gaussian form in the above equation, and several of 
its variant options, would be widely used in the years to 
come, mainly by the experimentalists, due to the simplic- 
ity of the emission source and analytical results allowed 
by this profile. We will see which are the parameters and 
interpretations derived from it in a while. 

1.3 SIMPLE PICTURE 




FIG. 3: Simplified picture: two point sources, I and II, emit 
quanta considered as plane waves, which are observed in de- 
tectors A and B, respectively, with momenta k'^ and kt^ ■ Since 
the quanta are indistinguishable, there are two possible combi- 
nations for this observation, illustrated by the two continuous 
and the two dashed lines. 

The amplitude for the process can be written as 
v2 



where the (-f) sign refers to bosons and the (— ) one, to 
fermions. In the above equation, (pi corresponds to an 
aleatory phase associated to each independent emission 
(completely chaotic sources), i.e., one phase at random 
in each emission. These phases are also considered to be 
independent on the momenta k of the emitted quanta. 

The probability for a joint observation of the two 
quanta with momenta ki and k2 is given by 

P2(A:i,A:2)- (|A(fci,fc2)|') = 

= -[2± (e*('^'i"'''2)-(2:i-a:2)^g±i(0l-H02-0i-02)^ -l-C.C.)] 

= 1 ± cos[(fci - fc2).(a;i - 2:2)]. (3) 

The emission being chaotic, we have to consider an 
average over random phases, i.e.. 



At this point, it is natural to ask the question: How 
to understand interferometry, or two-particle correla- 
tion, in a simple way? First of all, we should antici- 
pate that it follows from considering two essential points: 
the adequate quantum statistics and chaotically emitting 
sources, which was already emphasized by Bartknik and 
Rz§,iewski|3- Let me illustrate it by a simple example 
of only two point sources, as shown in Fig. 3: 



The two-particle correlation function can be written as 

C(fci,fc2) = -^2(^1,^2) ^ l±cos[(fci-fc2).(a;i-2:2)] , 
-Pi(ki)Pi(K2) 

(5) 
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where Pi{ki) is the single-mclusive distribution. It is es- 
timated in a similar way as in the simultaneous detection 
discussed above, i.e., 



-iki .(x A-xi) ^i4>i ^-iki.{xA-X2) 



(6) 

In the above case, we would have (e='=«(0i-02)^ _ <50^02- 
Since the source is supposed to be chaotic, the two 
aleatory phases of emission would be equal only if they 
were emitted at the same space-time point. However, 
since we are considering here that the probability of two 
simultaneous emission by the same source is negligible, 
we would be forced to conclude that only possible solu- 
tion to this problem that would satisfy this criterium is 
that the average over phases is null, in the case of obser- 
vation by a single detector. We see then that Pi{ki) = 1 
in this case and then the result on Eq.jS)) follows. 

Already from the very simple example discussed above, 
se can see that, in the case of two identical bosons 
(fermions), we expect to see that C{q = fci — ^2 = 0) = 
2 (0) for completely chaotic sources. On the contrary, in 
the case of total coherence C{q = fci — ^2) = 1 for all 
values of the momentum difference. For large values of 
their relative momenta, however, the correlation function 
should tend to one, which is clearly not the case in Eq. 

. But this is merely the consequence of considering an 
oversimplified example of only two point sources. 



1.4 EXTENDED SOURCES 

More generally, for extended sources in space and time, 
if p{x) is the normalized space-time distribution, we have 

= Pl(fcl)Pi(fc2) J d^Xi J d^X2 \A{ki,k2)fp{xi)p{x2) 

= Pi(fci)Pi(fc2)[l±|p(g)P], (7) 
where 

p{q) = [ d^x e^^'^^pix) (8) 



is the Fourier transform of p{x). Conventionally, we de- 
note the 4-momentum difference of the pair by q'^ = 
(fc^ — k!^), and its average by K'^ — ^(fc^ + kl^). 

Then, the two-particle correlation function can be writ- 
ten as 

^(^1' = p?f^ p^rl^ ^ - 1 ± ^ 1/5(9)1'- (9) 
Pi(ki)Pi[k2) 

In Eq.© we added, as historically done, the parameter 
A, later called incoherence or chaoticity parameter. 



This was introduced by Deutschmann et al.[3|, in 1978, 
as a means for reducing systematic errors in the experi- 
mental fits of the correlation function. The origin of the 
large systematic errors was the Gaussian fit. The reason 
was that the experimentalists tried to fit the data points 
with Gaussian functions whose maxima in q = were 2, 
although the data never reached that maximum value. 
This led to discrepancies and to large systematic errors. 
The easiest way out of this apparent inconsistency was 
to add a fit parameter, A, thus reducing the systematic 
errors by the introduction of this extra degree of freedom. 
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FIG. 4: Simple illustration corresponding to the ideal Gaus- 
sian source. The upper curve represents to the bosonic case, 
while the lower curve, the fermionic one. The parameter R is 
the r.m.s. radius of the emitting region. 

To illustrate the correlation function as written in 
Eq.® with a simple analytical example, let us consider 
the Gaussian profile, i.e., 



p{x) 



P(9) 



(10) 



Consequently, in this very simple example, a typical cor- 
relation function is written as 



C(fci,fc2) = l±Ae 



(11) 



In equation l(TT|l. as we denoted before, the plus sign 
refers to the bosonic case, and the minus sign to the 
fermionic one. We easily see that, in this simple example, 
we would expect experimental ideal HBT data to behave 
as sketched in Fig. 4, where the upper part refers to 
bosons and the lower one, to fermions. We see that, in the 
two-boson (two-fermion) case, there is an enhancement 
(depletion) of the correlation function in the region where 
the relative momenta of the pair are small. In both cases 
of this simple example, the typical size of the emission 
region corresponds to the inverse width of the C{k\,k2) 
curve, plotted as a function oi q = k\ — k2. 
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Returning to the discussion of the fit parameter A, f 
would hke to point out that there is a very simple ex- 
planation to reconcile this apparent inconsistency, with- 
out the need to introduce this extra degree of freedom. 
Limited statistics is behind it, since it is virtually impos- 
sible to measure two identical particles with exactly the 
same momenta. This led the experimentalists to split the 
momenta of the particles in small bins. In more recent 
times, these bins can be projected in two or more dimen- 
sions. For instance, along the income beam direction in 
fixed target heavy ion collisions (ql), and in the direction 
transverse to it (qt)- Good quality data allow the experi- 
mentalists to consider very small bin sizes. Nevertheless, 
their range is finite. Being so, when the correlation func- 
tion is projected along, say, the qr direction, the small- 
est value of ql is not zero, but within the first (smaller) 
bin size, in case of high enough statistics. Consequently, 
we immediately see that the correlation function plotted 
as a function of qr, will not reach the maximum (min- 
imum) value of 2 (0) for bosons (fermions) at qr = 0. 
Naturally, the larger the first bin size, the bigger is the 
deviation from the maximum (minimum) expected value 
for the correlation function at g = 0. This can be bet- 
ter seen with the help of Fig. 5, where the two upper 
curves represent the pure correlation function and the 
two lower ones, the theoretical correlation functions cor- 
rected by the Gamow factor, T(q), a multiplicative factor 
taking into account 2-body Coulomb final state interac- 
tion. This factor distorts the pattern, mainly at small 
values of the momentum difference of the pair, and is 
written as 



T(<?) 



qc/q 



exp(gc/g) - 1 



= 27ram 



(12) 



where q is the four-momentum difference and a is the fine 
structure constant. The Gamow factor simply multiplies 
the entire expression in Eq.©, Hll|l . and all other forms 
of correlation function for two charged, identical parti- 
cles. Fig. 5 was generated by the code CERES, whose 
hypothesis and formulation will be discussed later in this 
manuscript, in Section II. 2. 

We should emphasize, however, that the simple rela- 
tion between the two-particle correlation function and 
the Fourier transform of the space-time distribution, as 
written in Eq. Hll|l is not straightforward, in general. It 
is observed only when we can consider the phase-space 
as decoupled, i.e., /(x, p) — p{x.)g{p), where p(x) is the 
source space-time distribution, and ^(p) is the energy- 
momentum distribution. In general, it cannot be decou- 
pled in this way. This is, for example, the case in rel- 
ativistic heavy ion collisions, where the source expands 
during its life-time. The reason is that HBT is not only 
sensitive to the source geometry at particle emission but 
is also sensitive to the underlying dynamics. This makes 
the analysis model-dependent and more powerful for- 
malisms (like the one proposed by Wigner, the Covariant 
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FIG. 5: The plot illustrates the fact that we can have the 
maximum of the bosonic correlation curve below unity, even 
when the intercept parameter is fixed to be X = 1. This nor- 
mally happens in realistic cases of finite statistics. In the plot, 
we see the two-pion correlation as a function of qr, for two 
possible bin sizes, i.e., for qL < 0.01 GeVJc and qL < 0.1 
GeV/c. This plot was extracted form Ref.]jj 



Current Ensemble, etc.) must be adopted. We will dis- 
cuss more about this limitation later on. 



1.5 FURTHER APPLICATIONS 

In the 1970's, Kopylov, Podgoretskii, and GrishinQ 
used second-order interferometry to study several inter- 
esting problems. For example, they modelled the nucleus 
as a static sphere with radius R, emitting pions from its 
surface and got the following correlation function 



C(fci,fc2) 

where 

qx 

9o 



(13) 



= q 



K 



K| 

q-qii 



El — E2 

1 
2m 



1 



-(k?-k^) 



2m 

ki - k2).(ki + k2) cx q\\ 



(14) 



The two variables in Ea. (|14|l are nowadays known 
as Kopylov variables. With respect to the same 
parametrization, CocconiQ re-interpreted in 1974 the 
quantity ct = 5{< R) as the thickness of the pion emis- 
sion layer. They used similar forms for studying: i) the 
lifetime of excited nuclei through the interferometry of 
evaporated neutrons, ii) shape and size of multiple pro- 
duction region with 1:^1:^ correlations, etc. They also 
applied to CERN/ISR data on pp,pp. 
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Many other scientists contributed to the field during 
that decade. Just to mention a few names, I would 
quote Shuryak; Biswas; Fowler & Weiner; Giovannini 
& Veneziano; Grassberger; Yano & Koonin; Gyulassy, 
Kauffmann & Wilson, etc. Many of these contributions 
were organized in a collection of reprints, edited by R. M. 
Weiner |10|. in the late 1990's, which is a very good source 
of these reference papers. Among them, the papers by 
Gyulassy, Kauffmann & Wilson jllj, as well as those by 
Fowler and Weiner 12], represented important steps in 
the field, for introducing more powerful formalisms for 
studying the cases of coherent, chaotic, and partially 
coherent sources. On the other hand, Grassberger [l3j 
called the attention to the fact that resonances could 
play an important role in interferometry, since the long 
lived ones could distort the correlation function in the 
region where it was more significant, i.e., at small val- 
ues of q, thus changing the chaoticity parameter con- 
siderably. This is due to the fact that a resonance of 
4-momentum k, mass M, and width F, would travel a 
distance ~ k/{MT) before decaying, causing interference 
effects whenever k.q < MT. The first attempt to analyze 
the effect of resonances on interferometry in detail was 
made more than a decade afterwards. I will discuss it 
later in Sec. n.2. 

Despite the comment made earlier in the text, regard- 
ing the limitations of the static Gaussian fit, this has 
always been the preferred source model, due to its sim- 
plicity. Along this line, it is instructive to observe that 
there is a Gaussian limit of Ea. p3|l . corresponding to 



, (15) 



a variation of which suggests a non-relativistic 
parametrization of the correlation function, i.e., 

C(fci,fc2) - 1 ± A e-«o-'/2e-«-«-/2e-«'^-/2 . (16) 

The expression in Eq. (|16() has been widely employed 
since the beginning of high energy heavy ion collisions, 
becoming the standard form to analyze two-particle in- 
terferometry, particularly among the experimentalists. 
In Eg. 1)16(1 . ql is the momentum difference along the di- 
rection of the incident beam, qx is the component trans- 
verse to beam direction, and is the time component. 
In the late 80 's, the most popular form changed slightly, 
according to the suggestion by Bertsch0|, becoming 

C(fci, fca) = 1 ± A e-«s«le-9o«?.e-«' /2 ^ (^7) 

similarly to the previous definition. However, in 
Bcrtsch's suggestion, there was a decomposition of the 
transverse component, partially incorporating the defi- 
nition introduced by Kopylov and Podgoretskii, i.e., qo 
and qs are both perpendicular to beam direction but 



qo II Kt [= ^(kiT +k2T)], and qs -L Kx- As before, qL 
represents the component of the pair momentum differ- 
ence along the beam direction. Latter, Heinz et al., sug- 
gested to include a out-longitudinal cross term in Gaus- 
sian fits to the data, i.e., the correlation in this case would 
be written as[l5| 

C{kuk2) ^l±\e-'^"sRle~<'lRle-ilRl^-2ql,RhL . 

(18) 

Ever since, this field has been under constant develop- 
ment and expansion, both in the theoretical and in the 
experimental grounds. I will briefly highlight only some 
of the theoretical contributions to the field, mainly fo- 
cusing at the ones from the Brazilian group and some 
collaborators from abroad, since a complete discussion 
of the contribution along theses 50 years is beyond the 
scope of the present review. For a complete survey of 
the subject, as well as of the theoretical and experimen- 
tal progress in the field, I would stron g ly enco urag e the 
reader to look into Refs. O E El El 111 US . 



II. CONTRIBUTIONS FROM GROUP MEMBERS 

The first contact of the group, whose contributions we 
are discussing here, with HBT interferometry started in 
the mid- to late eighties, and was the subject of my PhD 
Thesis ll^. In fact, it happened a few years before the 
group itself began performing as a group. Nevertheless, 
this topic consistently appeared during the group meet- 
ings along these years and, since this is a historical per- 
spective, it is worthy to insert the subject in this context. 
Around the beginning of the decade of 1980, there was 
already an emerging subject that was attracting the at- 
tention of the high energy community: the possible exis- 
tence of a new state of matter, the Quark-Gluon Plasma 
(QGP), expected to be produced at high enough tem- 
peratures and/or densities. The QGP is a state in which 
quarks and gluons, the constituents of the hadrons, would 
be free to wander around a volume much bigger than the 
usual hadronic size. This state was expected to exist for 
a brief period of time, since only usual hadrons, with 
quarks and gluons confined in their interior, have been 
observed empirically. This imposed the need to look for 
probes of its existence. Among them, Interferometry was 
suggested, as a means to estimate the dimensions of the 
system formed in high energy collisions, thus testing if 
it was produced in such a new state of matter. In fact, 
James D. Bjorken was the person who suggested pion 
interferometry as the subject of my Ph.D. thesis |2^. 



II.l EXPANSION EFFECTS IN HBT 

In the first paper on the subject, we started by mak- 
ing the hypothesis that the Quark-Gluon Plasma was 
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already being produced in pp and pp collisions at the 
CERN/ISR. We considered |25l| that the system produced 
in such collisions expanded before emitting the final par- 
ticles (hadrons), according to the one-dimensional Lan- 
dau Hydrodynamical Model l^^. In the initial stage, 
the system was formed in the QGP phase at a certain 
temperature, Tq, started expanding and cooling down, 
until it reached the critical temperature, Tc, which we 
assumed to be of order of pion mass. It could be imag- 
ined that, once Tc was reached, the hadronization oc- 
curred instantaneously, followed by the particle emission. 
This simplifying hypothesis was actually adopted in the 
general study of the effects on the correlation function 
caused by the system expansion. On the other hand, 
the energy density of the ideal QGP fluid once Tc is 
reached, is much higher that the correspondent one for 
a hadronic system, due to the statistical degeneracy fac- 
tors. More explicitly, eqapiTc) = ^{gg + lgq)T^ + B, 
where gg = 8{color) x 2{spin), gq — 2{qq) x 'i[color) x 
2(spin) X Nf (flavors) are the gluon and quark degen- 
eracy factors. The constant B is the vacuum pres- 
sure in the MIT Bag model. And, the hadronic cor- 
respondent for an ideal gas of pions and kaons, at Tc, 
is e^(Tc) = ^[{g-^)(l){m^/Tc) + gK'Pi'mK/Tc)], being 
3 and gx = 4, respectively the statistical factor 
for pions and for kaons (The function ^(z) is a combina- 
tion of Bessel modified functions of second class, Ki{z), 



quarks and giuons, with initial radius R R 



proton • 
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T.::?=o(±ir{ 



3K2[z(l+m)] 
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Nevertheless, the large ratio of the QGP to the hadronic 
statistical degeneracy factors, together with the entropy 
conservation during the phase transition, make the dura- 
tion of the mixed phase very long. And mesons would be 
emitted during all that period. This was a more realistic 
hypothesis that was adopted when comparing our pre- 
dictions with experimental data. However, for the sake 
of simplicity, we considered in the calculation that the 
emission occurred at a typical average freeze-out time, 

< T/ >. 

In our calculation, we neglected the transverse expan- 
sion and used the asymptotic Khalatnikov solution, i.e. 



a 



To 



2 \t-x 



In 



(19) 



where a is the system rapidity, cq ~ 1/VS is the sound 
velocity \/(l — Cq)/ttI, being / the initial thickness of the 
fireball (solution valid whenever ^ ^ |a|). This is essen- 
tially the Bjorken picture of hydrodynamics but with dif- 
ferent initial conditions. In the version we adopted of the 
hydrodynamics, the initial temperature Tq depends on 
the value of the fireball, with mass M, initially formed in 
high energy hadronic collisions. The hypothesis we used 
was that a large Lorentz-contracted fireball was formed 
around one of the incident particles 5?), constituted of 



fireball mass was estimated as the missing mass, i.e., by 
discounting the fraction of the energy available in the 
center of mass of the p{p)p collision that was dragged by 
the dominant particle after the collision happened. For 
relating this initial temperature with the fireball mass, 
we equated the number of produced hadrons at Tc to 
the (conserved) entropy, which can easily be estimated 
for a QGP as S{To) = s{To)Vqgp. The initial volume 
of the QGP can be related to the initial energy density, 
eq, which can be simply written as eq = M/Vqgp- The 
initial entropy density, s{To), can be estimated through 
statistical relations, leaving the calculation of the initial 
QGP volume to be made. Since we considered p{p)p colli- 
sions, we assumed that Vqgp = |7ri?Q^^, where Rq and 
TOp are, respectively, the proton radius and mass. At 
the end, the final proportionality coefhcients were esti- 
mated with the help of the experimental data on charged 
multiplicity versus the missing mass. Finally, the initial 
temperature was related to the fireball mass by a numer- 
ical factor. To = 0.0989-\/M. From that, we estimated 
Tc, instant corresponding to the beginning of the phase 
transition, Tc = A i and also the instant it ended, 

Tf, as well as the typical (average) duration, < Tf > (see 
Ref. m for details). 

We had adopted the Kopylov variables described be- 
fore, in Eq. H14|l and sketched in Fig. 6, as the relevant 
momentum difference of the pair of pions. 




FIG. 6: Illustration of the Kopylov variables: we see that 
qi = \\ K and qr = q1 -L K . The figure also shows other 
notations commonly used: q = Ap and K = p. 

For studying the general behavior of the correlation 
function under the influence of the expansion effects, we 
assumed that each point on the surface r = r/ of the 
QGP, where T = Tc ~ m^, was an independent chaotic 
source with momentum spectrum given by 



f{p) 



1 



(27r) 



E 



exp (- 



Tc ' 



(20) 
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where u'^ = (cosha, sinha, 0, 0) is the 4- velocity of the 
fluid, and = {E,px,Py,Pz) is the 4- momentum of the 
emitted particle. The amplitude for a particle emitted at 
x' to be observed at x is written as 

A{x,x') = f dp^/mexp[-ip^{x^ - x'nV'^^^"' , 



(21) 

where (f>{x') is a random phase. We followed the formula- 
tion and notation of Ref . |28l | for writing the probability of 
detecting two quanta of momenta pi and p2 in an event, 
as 

W{pi,P2) = /(0,pi)/(0,P2) + -P2), \iPl +P2)]\^ , 

(22) 

where 

l[(Ap),.|./...(A.,»*— ./..'.,,A.,.',, 

(23) 

and 

At At 
(A*[x- ^,x']A*[x+^,x']) ^Six' -x)I{x,Ax,x') . 

(24) 

The average indicated in the above equation is taken over 
the random phases 4>{x') and 4){x"). In particular, we see 
that the single-inclusive distribution is written as 



W{pi) = (I / dx' / dxe'^^-'''A(x,x')dx\'^) = /[0,p, 



(25) 



' I I I ' I ' I M ' 

M = 540GeV 
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From this formulation we obtained the expression for 
the two-particle correlation function. With that, we stud- 
ied several different kinematical zones, trying to appre- 
hend the lessons that idealized theoretical cuts could 
teach us. Many interesting and important results came 
out of that study (see Ref. [i^ for details). 

The first, unexpected effect, was showing a clear influ- 
ence of different emission times at which the pions were 
emitted from the source. I will write it in a simple analyt- 
ical form, later in the text. This effect showed itself into 
the correlation function when plotted as a function of q^, 
the Kopylov variable parallel to the average momentum 
of the pair, K = ^(ki + k2). It is illustrated in Fig. 
7. Independently of our knowledge, S. Pratt had also 
suggested that the time would influence the correlation 
function, so that large short-lived sources could result 
into a similar correlation function as a short long-lived 
one [igf . 

The dependence on the average momentum of the pair, 
Kt, shown above, was a symptom of effects coming from 
the underlying dynamics, and reflected the break-up of 
the naive picture, in which the correlation function de- 
pended exclusively on the variable qx, as in the Gaussian 
example discussed above. 
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FIG. 8: In (a), the two-pion correlation function is shown in 
as a function of the Kopylov variable qr, averaged over the 
interval ql < 0.15 GeV (we adopted natural units, in which 
h = c = 1). In (b), it is shown as a function of the other 
component, qL JJ K, for qr < 0.15 GeV. The data points 
are from Ref. fS^ J. These plots were originally published in 
Ref^. 



FIG. 7: C{pi,p2) as a function of the momentum difference 
ApT (= qT , in the text), for several values of the average mo- 
mentum p± (= Kt in the text). A clear dependence is seen: 
the curves broaden with increasing average momentum, show- 
ing a progressively smaller effective emission region. This plot 
is a reproduction of the one originally published in Ref. '25]. 



We also compared our results with data of pp and pp 
collisions at CERN/ISR (^s = 53 GeV) and could ef- 
ficiently describe the trend of data. This was maybe 
the only successful description of that particular experi- 
mental result, reflecting the need for more powerful for- 
malisms when describing HBT interferometry at high en- 
ergies. In this case, in an effort to make the estimate more 
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realistic, we considered that the emission occurred later, 
at a typical instant of time < Tf >, averaged over the long 
period that lasted the first order phase transition. We 
should notice that the curves in Fig. 8 are not fits, but 
predictions from the model, obtained without the need 
for introducing the A parameter (which is equivalent to 
fixing A = 1). Similar treatment within the same model 
was also given to two-kaon interferometry data from the 
same experiment |2ItI| . successfully describing data. 

Fig. 7 and 8 above also show another important result 
from that study: the observation of strong distortions 
in the correlation function, definitely departing from the 
Gaussian shape, due to the dynamical effects related to 
the expansion of the system. 
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II. 2 NON-IDEAL EFFECTS 

More than 10 years after Grassberger pointed out the 
important role resonances could have in interferometry 
we investigated it in detail, in collaboration with M. 
Gyulassy 30J. In particular, we analyzed the effect of 
resonances decaying into pions, following the predictions 
of resonance fractions from the ATTILA version of the 
Lund modeljsj. Very briefly, it can be understood as fol- 
lows: long lived resonances, such as w, 77, 77', can mimic 
sources with longer life-times, even if they freeze-out si- 
multaneously as the direct tt's. 

I call the attention to the fact that, although denoted 
by the same letters, the variables ql and qr of momen- 
tum difference, appearing in Fig. 9, are defined as the 
components parallel and perpendicular to the incident 
beam direction, respectively, as became a convention in 
high energy heavy ion collisions. 

It was very important to verify if we could explain the 
preliminary results obtained by the NA35 collaboration 
(CERN/SPS)f^ colliding 0+Au at 200 GeV/nuclear in 
a more conventional way, without assuming QGP for- 
mation. This last hypothesis had been suggested by 
G. Bertschjs^ at that time as the only possible expla- 
nation. Actually, by introducing resonances decaying 
into pions explicitly into interferometry, we managed to 
demonstrate that our result, assuming the formation of 
a regular hadronic system and no QGP, led to an equally 
good agreement with data, as can be seen from Fig. 9. 
In fact, we showed in Ref. 1301 (and presented in the 
Quark matter '88 Conference |3J|) that the NA35 pre- 
liminary data of that time was consistent with a wider 
range of pion source parameters when additional non- 
ideal dynamical and geometrical degrees of freedom were 
incorporated into the anal ysis by extending the Covari- 
ant Current Formalism [ill hd, l35l. 
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FIG. 9: The two-pion correlation function is shown in three 
different scenarios, as a function of qr ■ The first, in parts 
(a) and (d) were calculated assuming an ideal inside- outside 
cascade source (IOC) as in the 1-D Bjorken victure\3(\I. In 
parts (b) and (e), a non-ideal resonance gas source is consid- 
ered with parameters suggested by the ATTILA version of the 
Lund model. Parts (c) and (f) correspond to the QGP model 
of Ref.JSSiJ. Plots in the upper panel were calculated in the 
central rapidity region (^2 < < 3j, whereas the ones in the 
lower panel refer to 1 < < 2. This figure is a reproduction 
of the one originally published in R,ef.J[S(\l. 



II.2.1 Ideal Bjorken IOC Picture 

In the Covariant Current Ensemble formalism^, 13(1 

I35I1, the correlation function for identical bosons (since 
mainly bosonic HBT will be discussed in this review) 
can be expressed as 

C(fci,fc2)^g(g,j^)^l + A '^^'^'^f , (26) 

G{ki,ki)G(k2,k2) 

where qt" ^ k^ - fc^ and Rf" = Kfc^* + fc^). 

The complex amplitude, G(fci, fc2), can be written as 



G(fci,fc2) = J d^xd^pe"^ ''-Dix,p)foiu'}kiMu^fk2^) , 

(27) 

where D(x,p) is the break-up phase-space 

distribution [3(1 l37j and the currents, jo{uf.ki), contain 
information about the production dynamics. The 
one-particle spectrum is obtained from Eq. H27(l by 
imposing fci ~ k2, which leads to 



(28) 



G{h,h) = d xd p D{x,p)\jo{u^fki^)Y 
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The currents jo{uf.ki) m Eqs. (|27I28|I can be associ- 
ated to thermal models, and can be written in a covariant 
way as 



jo(fc) oc v/u^fc^expl ^} 



(29) 



However, to make the computation easier, we adopted a 
more convenient parametrization 



joiu.k) = exp{--— ^} , 



(30) 



where the so-called pseudo-temperature T,,., is related 
with the true temperature T according to 35] 

Tpsix) = 1.42r(a;) - 12.7 MeV . (31) 

This mapping between T{x) and Tps{x) was later shown 
to be a good approximation also in the case of kaon 
interferometry |4J| . 

With the covariant pseudo-thermal parametrization as 
in Eq. 1)30(1 . the complex amplitude can be rewritten in 
a simpler form. 



0(^1,^:2) = (e^'^-^e-^P-^/f'"^)) , 



(32) 



where the bracket < ... > denote an average over the 
pion freeze-out phase-space coordinates. 

In the case of Bjorken ideal Inside-Outside Cascade 
(IOC) picture, the phase-space distribution involves a 
fixed freeze-out proper time rj and a perfect correlation 
between 77 and y. The correspondent phase-space distri- 
bution is written as 



11 e 
D{x,p) = ■TrP—Siy-'r)Sirj-y)dipo-Ep)g{pT) 



where Ep = 



(33) 

1? is the energy and 5(pt) is the 
transverse momentum distribution; the rapidity distri- 
bution is considered to be uniform, i.e., ^ = p. In the 
ideal IOC picture, there is a perfect correlation in phase- 
space between the space-time rapidity 



and the energy-momentum rapidity, 

y = \ > 

2 E - pz 



(34) 



(35) 



i.e., they are indistinguishable. 

To obtain simple analytical equations, we assume a 
very narrow distribution of pt around small momenta, 
i.e., gipr) — 5'^{pt)- The finite pion wave-packets gener- 
ate the finite px distribution in our case. By substituting 
D{x,p) from into (|77l) and considering the pseudo- 
thermal parametrization (|30|l for the currents, the func- 
tion G{ki,k2) was found to be|35l| 



G{kuk2)^2<^>{^—,h{qTRT)}KoiO , (36) 
dy qrUT 

where 

= [-T^{miT -\- m2T) - 'iT{miT - m2T)f + 

ps 

2(7^+^'W"^2T[cosh(Ay)-l] , (37) 

ps 

and Ay = yi- y2. 

The single-inclusive distribution is then written as 



, N d^N dN ,171. 

Gih, h) = E^ ^2< — >Ko{ 



dyi 



-) 



(38) 



ps 



To compare theoretical correlation functions with data 
projected onto two of the six dimensions, we computed 
the projected correlation function trying to mimic what 
is done in the experiment, i.e.. 



CprojiQT, Ql) 



J d^kid^k2P2iki,k2)A2iqT,qL;ki,k2) 



Jd3k,d3k2Pl{k,)P,{k2)A2{qT, Ql; h, k2) 

(39) 

where Pi and P2 are, respectively, the single- and two- 
pion inclusive distributions. It is essential to have per- 
fect correspondence between the experimental informa- 
tion and the theoretical estimates concerning the num- 
ber of dimensions into which the correlation data is pro- 
jected, the cuts in momentum and rapidity, the sizes of 
the bins, etc. A2 is the experimental two-particle binning 
and acceptance function, through which we approximate 
the theoretical estimate to the empirical cuts. 

We should notice that, experimentally, the two-particle 
correlation function in high energy collisions is obtained 
by measuring the following ratio 



Cexpi'^-) j) 



^ 3Xp 



ACexpihj) — Cexpih j)] 





2 

+ 


[AS(*,j)l 









(40) 



The numerator A{i, j)±AA{i, j) corresponds to the com- 
bination of pairs of identical particles from the same 
event, and the denominator, B{i,j) ± AB{i, j), repre- 
sents the background; Afexp is an experimental normaliza- 
tion. Historically, the background for identically charged 
pions have been the combination of unequally charged 
ones. However, later it was realized that 7r"'"7r~ could fre- 
quently come from the decay of resonances, which would 
distort the background and cause strange pattern for the 
correlation function built in this way. They soon realized 
that a better way to construct the background was to 
combine identically charged particles, but from different 
events. Another possibility used sometimes is a Monte 
Carlo simulated background, taking into account the ex- 
perimental cuts and acceptance. 
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II.2.2 Non-ideal IOC 

The non-ideal picture mentioned before referred to the 
underlying effects that would be important to incorpo- 
rated into the interferometric analysis, even restricting 
the attention to completely chaotic sources. For instance, 
the TT^ rapidity distribution at 200 AGeV was clearly 
not uniform, as assumed in the asymptotic Bjorken pic- 
ture, but would be better described by a Gaussian with 
width Yc « l-4j23- On the other hand, a large frac- 
tion of the pions could arise from the decay of long lived 
resonances, such as oj, K* , rj, etc, as was suggested by 
Grassberger^^. In coordinate space, the finite nuclear 
thickness, together with resonance effects, could lead to a 
large spread (Ar) of the freeze-out proper times, and to a 
wide distribution of transverse decoupling radii (Rt)- In 
phase-space, there is not the perfect correlation between 
the space-time and energy-momentum rapidity variables 
present in the ideal Bjorken picture. Instead, as sug- 
gested by the ATTILA version of the Lund model, they 
would be better related by a Gaussian with finite width 
A77. Besides, other correlations may have to be consid- 
ered if collective hydrodynamic fiow occurs, for instance, 
between transverse coordinates {x±) and transverse mo- 
mentum component {p±_). All these effects together were 
generically called as the non-ideal picture, which is equiv- 
alent to considering a more realistic picture than the one 
idealized by the Bjorken in his version of the 1-D hy- 
drodynamics. The phase-space distribution representing 
these effects together can be obtained from Eq. H33|l by 
replacing 



p — 6{Tf 



t)<5(?/ - y) 



At- 



■ exp 



exp 



2y2 



■ exp 




/2ttAt] 

Besides the modification in Ea. H41() . there is a major cor- 
rection to be added, i.e., the effect of long-lived reso- 
nances decaying into pions. This can be included in the 
semiclassical approximation jsol l34| . The pion freeze-out 
coordinates, a;^, can be related to the parent resonance 
production coordinates, a;^, through 



(42) 



where is the resonance four velocity and r is the 
proper time of its decay. Summing over resonances r of 
widths r^, and averaging over their decay proper times, 
we obtain, instead of Ea. H32|l the final expression [sol Is^ 

exp(i(7Xr — KUr/Tr) , 



/r- 



(1 — iqUr/Tr) 



(43) 



where /(7r-/r) is the fraction of the observed 7r~'s arising 
from the decay of a resonance of type r, and the temper- 
ature Tj. characterizes the decay distribution of that reso- 
nance. According to the Lund model, the main resonance 



contributing to the negative pion yield at CERN/SPS en- 
ergies are .f(^-/^) = 0.16, f{^-/K') = /[tt- /(,,+,,')] = 0.09, 
/(7r-/p) = 0.40, /(7r~/djrect) = 0.19. Although we in- 
cluded direct pions and the ones coming from p decay 
independently, they are hardly distinguishable, since p's 
decay very fast. 

All these effects combined were simulated in a Monte 
Carlo code, named CERES, which was also able to in- 
clude simplified subroutines which mimicked the experi- 
mental cuts and acceptance [ifl 13^ . Is^ . 
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FIG. 10; Numerical results for the pion correlation function 
versus Ay are shown in the central rapidity region, when the 
non-ideal effects are introduced, one by one. Thi s pl ot is a 
reproduction of the one originally published in Ref. pi/ . 

For explicitly demonstrating how the above results 
would show themselves into the correlation function, we 
derived in Ref. [s^l several analytical relations where 
the non-ideal effects were progressively introduced. We 
started with the ideal Bjorken 1-D picture. Then, instead 
of the instant freeze-out assumed there, we considered a 
spread (At) and, subsequently, all the other effects men- 
tioned above, one by one. For fixed qr = 0, the results 
as a function of the rapidity difference Ay (« ^l/tot 
for Ay <C 1), can be seen in Fig. 10. For generating 
these plots, we estimated the correlation function using 
the code CERES and adopting the following values for 
the parameters: Rt ~ 4fm, r/ « 4fm/c, fixed the inter- 
cept parameter A = 1, At = 4 fm/c; Yc ~ 1.4, Arj « 0.8, 
plus the above fractions f^r/r from the Lund model, when 
adding resonances. 

In the curves shown in Fig. 9, nevertheless, we had 
three sets of parameters, corresponding to the different 
models compared there: in the Lund resonance case, be- 
sides the fractions /^/r, we used Yc ~ 1.4, A77 « 0.7, 
Rt « 3fm, tj « 3fm/c, y* = 2.5 and fixed A = 1; for 
mimicking the QGP model of Ref.js^, we used no reso- 
nance, Tf = 9.0 fm/c, Rt = 3.3 fm, A77 « 0.76, assuming 
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Yc oo and A = 1. 

The results shown previously in Fig. 9 clearly posed 
another problem to the interferometric probes of high en- 
ergy heavy ion collisions: several very distinct dynamical 
scenarios could lead to approximately the same final cor- 
relation function and similar experimental HBT results. 
One possibility to discriminate among very different sce- 
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between those two models, helping to separate long-lived 
scenarios from those where HBT results were generated 
by the effect of resonances. 

The important role of resonances stressed above led 
us to investigate how would be their influence for much 
smaller systems at higher energies. For this, we looked 
into pp and pp data from CERN/ISRjs^- For simplicity, 
in this calculation, we assumed the ideal Bjorken picture 
to describe the freeze-out distribution. Surprisingly, the 
result ^3 turned out to be neither compatible with the 
absence of resonance nor with the full resonance fractions 
predicted by the Lund fragmentation model mentioned 
before. Instead, the data seemed to be best described 
by the scenario with about half the resonance fractions 
predicted by the Lund model, as can be seen in Fig. 12. 
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FIG. 11: Comparison of two-pion, in parts (a) and (h), and 
two-kaon projected correlation functions, in parts (c) and (d) 
is considered as a function of the transverse momentum differ- 
ence, Qt ■ The plots are calculated in the central rapidity region 
and with ql < 0.1 GeV. Solid (dashed) lines indicate corre- 
lations without (with) Coulomb distortions. Parts (a) and 
(c) correspond to predictions based on the Lund model lSJjl . 
and parts (b) and (d), to the plasma hydrodynamical model 
of Ref.(3Sj. The pion data is from Ref.^S&J. Thts figure was 
extracted from Ref. Jf3^. 

narios, suggested in Ref. js^, was to explore distinct 
freeze-out geometries by comparing pion and kaon in- 
terferometry. This suggestion was motivated by the fact 
that an entirely different set of hadronic resonances decay 
into pions than into kaons. In the first case, according to 
the ATTILA version of the Lund fragmentation model, 
long lived resonances such as ui, rj and 77', contribute to 
the final pion yield, whereas, in the second case, half of 
the kaons are produced by direct string decay, and the 
other half by the decay of K*. On the other hand, in the 
QGP model considered previously in Ref.js^l for compar- 
ison, the freeze-out geometry of all hadrons was expected 
to be about the same. In the case of pions, we saw from 
Fig. 9 that both cases led to equally good results as com- 
pared to the experimental points. In the case of kaons, 
then, an entirely different behavior would be expected. 
Indeed, we see from Fig. 11a more significant difference 



FIG. 12: Negative pion correlation in pp and pp reactions at 
CERN/ISR energies are shown as a function of Qi„v on the 
left plot, and as function ofqx (with ql < 0.15 GeV/c) on the 
right plot. Dashed, solid and dot-dashed histograms indicate 
the calculated correlation functions with, respectively, no, half 
and full Lund resonance abundances. The data points are 
from Ref. J[S!ij on the left, and from Ref. ISAl .on the right 
hand side. Both plots were extracted from Ref.j^. 

Meanwhile, we also derived a general and powerful 
formulation ,4lj , based on the Wigner density formalism, 
which allows to treat complex system, by including an 
arbitrary phase-space distribution (i.e., the momentum 
distribution, ^(p), and the space-time one, p(x), could 
be entangled) and multi-particle correlations. This for- 
malism corresponds to a semiclassical generalization of 
the n-particle phase-space distribution, in which it is al- 
lowed for a Gaussian spread of the coordinates around 
the classical trajectories, in order to incorporate minimal 
effects due to the uncertainty principle. Also, correction 
terms due to pion cascading before freeze-out were de- 
rived using this semi-classical hadronic transport model. 
Such terms, however, can be neglected if the mean free 
path of pions is small compared to the source size, or 
if the momentum transfers are small compared to the 
pion momenta. The main result of that investigation can 
be summarized by the following formula for the Bose- 
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Einstein symmetrized n-pion invariant distribution 

TL 

exp [i{kj — krj.).Xj\ 

a- j=l 

SA{kj~K^,Pj)) , (44) 
with the smoothed delta function given by 



SA{k-k',p-K) 



exp{[p-i(fc + fc')]V(2Ap2)} 



exp I 



(27rAp2)-3/2 

^{k-kT] 



(45) 



The brackets < ... > denote an average over the 8n 
pion freeze-out coordinates {xi,pi, ...,Xn,Pn}, as ob- 
tained form the output of a specific transport model, such 
as a cascade or the Lund hadronization model'si'l. In 
the form written in Eq. (|44|l and H45|) , this formulation is 
ideally suited for Monte Carlo computation of pion inter- 
ference effects. The smoothed delta function results from 
the use of Gaussian wave packets with widths Aa; and Ap, 
which depend on details of the pion production mecha- 
nism. The sum runs over n! permutations a — (tri, (t„) 
of the indices; x, k,p, ... denote 4- vectors and all momenta 
are on-shell. This is a generalization of the Wigner type 
of formulation proposed in Ref. '28'| and used in the pio- 
neer work of HBT in the group, discussed in Section II. 1 
and in Ref. [2^ l25j . As a special limit, we found out 
that, for minimum Gaussian packets, i.e., AxAp — 1/2, 
having Ap ~ mTpg, where m is the particle mass and 
Tps is the pseudo temperature of Eq. (|31|) we recovered 
the interferometric relations derived within the Covariant 
Current Ensemble formalism|llj. which, for the sake of 
simplicity, was adopted as a first approach to the study 
of non-ideal effects on Interferometry described above. 
An equivalent alternative way of expressing Eq. (|44|l is 

m 



P„(ki,...,k„) « {CN..)Y,X{DA[q{],cj,),K{],a,)] 

a j=l 



(46) 



with 



DA{q,K) ^ d^x d^pe''"'D{x,p)5A{q,P~K), (47) 



where D{x,p) can be given by, for example, Ea. (|33|l . and 
CN,n = N\/{N — n)!, where n = 2 in case of two-pion 
correlations, and N is the multiplicity of the event. 

An interesting simple point explicitly demonstrated in 
Ref. is the dependence of the effective transverse 

radius on the average momentum of the pair, which was 
already shown in Fig. 7, as one of the results of Ref.[25|. 
and also suggested in [23 • This dependence on Kt ap- 
pears through the time dependence of the emission pro- 
cess. The demonstration was done by means of a simple 



Gaussian example, as in Eq. (|16|) . For better under- 
standing it, we should recall the definition of the average 
4- momentum of the pair, K'^ ^ and their difference, , 
i.e., 

q^ = {k^,-ki^)-K^^]^{k>^ + kl^). (48) 
From the above relations, it immediately follows that 



q^'K^ = g"A'o-q.K = ^ = 



_ q-K _ qt.Kt + qhKL 



Ek 



Ek 



(49) 

where Ek = VK^ + was written for the sake of sim- 
plicity. 

Propagating the above result into the Gaussian corre- 
lation function, we get 



R 
R 



2 



2i?|. + (Ar)2(K|/£;|) 
R:i + {Arf{Kl/El). 



(50) 



The results on Eq. (|50|l show that the time spread of 
the source freeze-out generally enhances the effective size 
measured by interferometry. 

I should remark that several contributions and invited 
talks presented in international conferences in the period 
are being omitted here, due to the lack of space. I would 
address to the Quark Matter Conference proceedings for 
that, as well as the proceedings of the RANP Conference 
and of Hadron Physics. 



II.3 DISCRIMINATING DIFFERENT 
DYNAMICAL SCENARIOS 

The coincidental agreement with data of two opposite 
scenarios, such as the resonance gas and the QGP dis- 
cussed before, in Fig. 9, stressed the necessity of finding 
other means to more clearly discriminate among differ- 
ent decoupling geometries. Although the comparison of 
kaon with pion interferometry was shown to be helpful, 
as seen in Fig. 11, it still lacked from more quantitative 
information. Then, how to disentangle different models 
in a more precise way? 

In order to answer this question, M. Gyulassy and my- 
self developed a method, in which a 2-D analysis was 
performed, comparing two-dimensional theoretical and 
experimental pion interferometry results|4^. For illus- 
trating the method we performed the calculations using 
the code CERES mentioned above, for two very distinct 
scenarios. The first one considered the effects of res- 
onances decaying into pions, including Lund resonance 
fractions. The other one ignored the contribution of 
resonances. The data points were kindly sent to us by 
Richard Morse, from the BNL/E802 Collaboration, and 
corresponded to Si + Au collisions at 14.6 GeV/c|4a|. as 
measured at the BNL/AGS. I will summarize the method 
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E802 Data No-res R=4,6,t = 3. 4 Full-res R=3,1 ,t=1 .6 
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FIG. 13; Tv~TT~ correlations in central Si+Au collisions is 
shown as a function of qr and ql. The preliminary E802 
data^'ASj corrected for acceptance and Coulomb effects are 
shown in part (a). Parts (h) and (c) show theoretical cor- 
relation functions filtered with the E802 acceptance. They 
correspond, respectively, to cases without and with resonance 
production. This figure was extracted from Ref.JJ^. 

by recalling Eq. ^ and Eq. The E802 experi- 

mental acceptance functions for two particles was approx- 
imated by 

A2iqT,qL;ki,k2) = Ai(fci)Ai(fc2)e(20 - |0i - 02|) 
SiqL - \k,i - k,2\)6iqT ~ I^Ti - fcrj) ■ (51) 

The angles are measured in degrees and the momenta in 
GeV/c. The single-inclusive distribution cuts were spec- 
ified by 

Ai(fc) = 6(14 < < 28) 

Q{viah < 2.2 GeV/c)e{y,mn > 1.5) . (52) 

The input temperature matching the experimentally ob- 
served pion spectrum was T « 170 MeV. 

For the purpose of performing a quantitative analysis 
of the compatibility of different scenarios with data, we 
computed the goodness of fit on a two-dimensional 
grid in the {qT, 9l) plane, binned with 6qT = Sq^ = 0.01 
GeV/c. The variable was computed (as sug gested by 
W. A. Zajc) through the following relation|42j 

^ ^''-^^ {[AA{^,JW + [J^^-'CtH{^,3)AB{^,JW} ' 

(53) 

where Af^ is a normalization factor estimated as to min- 
imize the average , which depends on the range in the 
9T)9l plane under analysis. The minimization of the 
average x^ was performed by exploring the parameter 



space of the transverse radius Rt and the time r, and 
computing the < >, averaging over a 30x30 grid in 
the relative region (Jt,9l < 0.3 GeV/c of relevant HBT 
signal. In the vicinity of the minimum we determined the 
parameters of the quadratic surface 

(x^(i?T,r)) =xLn + «(^T-i?T„)' + /3(r-ro)2 . (54) 

TABLE 1: 2D-x^ Analysis of Pion Decoupling 



Geometry 



x'{Rt,At) 


No Resonances 


LUND Resonances 


E802 Data Gamow Corrected 




2.1 


2.2 


Rqt 


4.6± 0.9 


3.1± 1.3 


Ato 


3.4 ± 0.7 


1.6± 1.0 


a 


0.027 


0.014 


P 


0.042 


0.023 



The quantitative differences could be seen in a 3-D plot 
of the correlation function, projected in terms of qT and 
qL^ shown in Fig. 13. The main results in the case of 
Gamow corrected data (i.e., where the HBT signal has 
recovered at small values of q by multiplying by the in- 
verse of the Gamow factor, T{q)~'^ , defined in Eq. H12|) '). 
are shown in Table 1. For a more complete discussion of 
the method, I would address the full article, in Ref. 4^. 
We see from the upper plots in Fig. 13 that the 2-D cor- 
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FIG. 14: The plot shows the number of standard deviations 
from unity, of the average < > P^^ degree of freedom, as 
a function of the range Qmax of the analysis in the {qT,qL) 
plane. This plot was extracted from Ref. tl3{l. 

relation function for the non-resonance case is clearly dif- 
ferent from the resonance scenario but only at very small 
values of the momentum difference qx and q^. Never- 
theless, if we look into the lower panel where we plotted 
the x^ distribution of the theoretical curves compared to 
the experimental points, we see that the distinction is 
blurred by the large fluctuations of data, mainly at the 
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edge of the acceptance. The most efficient measure of the 
goodness of the fit in this case was obtained by study- 
ing the variation of the average < > per degrees of 
freedom in the {qTilh) plane with respect to the unity, 
i.e., I < Xmin > this way, we found out that 

the resolving power of the distinction between different 
scenarios was magnified, as shown in Fig. 14. 

The method was further tested later, in a more chal- 
lenging situation|44l|. by comparing the interferometric 
results of K'^K'^ of two distinct scenarios, i.e., Lund 
predicted resonances (only K*'s and direct kaons con- 
tribute significantly to this particle yield) with the non- 
resonance case. This was done by Cristiane G. Roldao in 
her Master Dissertation, under my supervision. The data 
on interferometry from 81 + Au collisions at 14.6 

GeV/c was sent to us by Vince Cianciolo, from E859 Col- 
laboration (an upgrade of the previous E802 experiment 
at BNL/AGS). 



TABLE 2: 2D-x^ Analysis of Kaon Decoupling 
Geometry 



x'(i?T,Ar) 


No Resonance 


LUND Res. 




Uk,.,. = 1) 


(/k,,„ = /k/k. = 0.5) 


Optimized Rt and Ar 


(Xmin)30x30 


1.03 


1.02 


(Xmin) 10X10 


1.17 


1.30 


Rto 


2.19± 0.76 


1.95 ± 0.89 


At,, 


4.4± 2.0 


4.4± 2.6 


Q 


0.0410 


0.0299 




0.0058 


0.0034 



The acceptance function for the E859 experiment was 
approximated^^ by 

(gr , gL ; A:i , fc2 ) = (fci )yli (A;2 ) e (22 - 1 01 - 02 1 ) 
S{qL - \kzi ~ k^2\)S{qT - \kTi - krj) . (55) 

The angles were measured in degrees and the momenta 
in GeV/c. The single inclusive distribution cuts are spec- 
ified by 

Ai{k) = e{U<9iab<28)e{piab< 2.9 GeY/c) 

e(y™„> 0.75) . (56) 

In the case of kaons, the input temperature matching 
the experimentally observed kaon spectrum was T « 180 
MeV. 

It was expected to be harder to differentiate both sce- 
narios due to the lack of contribution from long-lived res- 
onances. It was found that they could still be separated, 
with data favoring the non-resonance scenario at the 14.6 
Si + Au collisions (BNL/AGS). As in the two- pion inter- 
ferometric analysis, the variation of the average < > 



FIG. 15: Zone in the (Rt,At) plane investigated, leading to 
the determination of the most probable region where the min- 
imum (x^); associated to (Rt,/^t), could he located. Parts 
(a) and (b) correspond to cases where the contribution of K* 
were ignored. Parts (c) and (d) were estimated including 
their contribution to the kaon yield. In parts (a) and (c) we 
fixed Ar = 0, and optimized only Rt. Figure extracted from 
Ref.^. 

per degrees of freedom in the (gr, qh) plane was the sig- 
nificant quantity to look at, as can be seen n Fig. 16. The 
main fit results found in this analysis are summarized in 
Table 2. For more details, see Ref.0- 
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FIG. 16: The number of standard deviations from unity of the 
average < X^ > P^''' degree of freedom is shown for increasing 
number of bins. In part (a), the transverse radius parameter, 
Rt, and the freeze-out spread, Ar, used for generating the 
theoretical plots, were optimized. In (b), only Rt was opti- 
mized, whereas Ar was fixed to be zero. Figure extracted from 
RefU^. 

The tests also ruled out the possibility of a zero decou- 
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pling proper-time conjectured by that experimental re- 
sults of AGS/E859 CoUaborationji^. This can be seen 
from part (b) of Fig. 16: by fixing Ar to be zero the 
numerical deviations of the average < > P^r degrees 
of freedom are completely meaningless. 



II.4 SONOLUMINESCENCE BUBBLE 

In the beginning of this review, we saw that the HBT 
interferometry was originally proposed for measuring the 
large sizes (of order R ~ 10^° m) of stelar sources in 
radio-astronomy . On the other hand, in the totality 
of the cases discussed here so far, the dimensions went 
down to the order of the hadronic or to the nuclear size 
(roughly, i? 10^^''' m). In between these two very dif- 
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FIG. 17: The upper part of the figure shows correlation func- 
tions plotted as function of Au) = loi — UJ2 ~ c(ki — ^2). 
Cases A to E correspond, respectively, to the first up to the 
fifth examples in Table 2. The lower part, shows the ge- 
ometrical form factor "l>(g) plotted as a function of X = 
■\/— (l/2)[d2$(0)/dg2]g. Figures extracted from Ref. 



ferent scales, Yogiro Kama, Takashi Kodama, and myself 
[4^ discussed in 1995 an interesting approach to a beauti- 
ful problem, in a distinct environment. The focus was in 
a small sonoluminescent bubble, whose radius would lie 



in the range R ^ 10^'^m. The phenomenon had been dis- 
covered long time ago, in 1934, at the Univ. of Cologne, 
but its single bubble version was found out by Gaitan et 
al.^3 O'^ly 1988. In this last case, a single bubble of 
gas (usually filled with air) formed in water, is trapped 
by standing acoustic waves, contracting every 10-12 pico- 
sec approximately, and simultaneously emitting light. In 
other words, the sonoluminescence process converts the 
acoustic energy in a fluid medium into a short light pulse 
emitted from the interior of a collapsing small cavitat- 
ing bubble [i^. The spectrum of emitted hght is very 
wide, extending from the visible to the ultraviolet re- 
gions. An important fact is that the light emission takes 
place within a very short period of time. The small size of 
the emission region and the short time scale of the emis- 
sion process make it difficult to obtain precise geometrical 
and dynamical information about the collapsing bubble. 
Besides, a particular aspect of the emission process is still 
controversial: some authors attribute the light emission 
to the quantum-electrodynamic vacuum property based 
on the dynamical Casimir eff ect l48l| , whereas others con- 
sider that thermal processes [43, such as a black-body 
type of radiation, should be the natural explanation. In 
Ref. [iil ■ we suggested that two-photon interferometry 
could shed some light to the problems raised above, by 
estimating the very small size and life-time of the sin- 
gle bubble sonoluminescence phenomenon. In fact, the 
simple existence of a HBT type of signal from such a 
collapsing bubble selects the scenario of explanation be- 
tween the two classes mentioned above. This is possible 
because the emission processes are opposite in the case 
of thermal models and in the Casimir based ones. In 
the first case, the emission is chaotic, a condition that 
is essential for observing the two-identical particle HBT 
correlations, and an interferometric pattern would be ob- 
servable. For the Casimir type of models, however, the 
emission is coherent and no HBT effect would be observ- 
able since, in that case, the correlation function acquires 
the trivial value C(fci, ^2) = 1 for all values of ki. 

In this work we neglect all the dynamical effects dis- 
cussed in the previous studies, and considered that the 
space and time dependence would be decoupled as well. 
We adopted the notation used in applied for our 
case. 



C{ki,k2) = 1 + 



1 



2 5(0,fci)5(0,fc2) 



(57) 



where 



S{q,K)^ / d^xe-''^-p{x)r{k,)j{k2) 



q — ki — k2 , K — (fci + k2)/2 , and j{ki) is the ampli- 
tude for the emission of a photon with the wave-vector 
ki at the source point x. The factor 1/2 in the second 
term of Eq. H57|l comes from the spin 1 character of the 
photon [sof. 
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We suggested a set of candidate models that could try 
and describe the system, as summarized in Table 3. 

TABLE 3 

Analytic expressions for some source parameterizations 
and the corresponding correlation functions are shown. 



Source p(r, t) 


Ciki,k2)-l 


A: e-'-'/2«^ e-*'/2-^ 




B: S{r - R) e-*'/^^' 


g-(A<.)-r- [sin{qR)/{qR)]^/2 


C: e{R-r) e-*'/2-' 


(9/2) e-^^^y^^^iqR)}-^ 
{ [cos{qR) ~ sm{qR)/{qR)]}^ 


D: e---/^ 9(3x2 - t2) 


[sm{ALuV3T)/ALuV3TY 

(l+g2i?2)-V2 


E: 

e{Rt-r) e-*'/^'9(t) 


9|/|V(8a*«), 

/ = + flZ+)W{z+)- 

{1- ^iz-)Wiz-)]-2fi 



As far as I know, no HBT experiment has been made 
yet with the photons emitted by the sonoluminescence 
bubble. T. Kodama and members of his group had initi- 
ated the experiment to produce single cavitating bubbles, 
with intension to perform the interferometric test I have 
just described above. In any case, the simple observation 
of a two-photon HBT correlation would be enough to rule 
out one of the two classes of models that aim at explain- 
ing the phenomenon, i.e., the ones based on the Casimir 
effect, in which the emission is coherent. It would auto- 
matically enforce the other class of thermal models, for 
which the emission is chaotic. 



II.5 CONTINUOUS EMISSION 

In Section II. 2.1, we discussed the Bjorken Inside- 
Outside Cascade (I0C)|3G| picture. It considers that, 
after high energy collisions, the system formed at the ini- 
tial time T — To, thermalizes with an initial temperature 
To, evolving afterwards according to the ideal 1-D hy- 
drodynamics, essentially the same as Landau's version 
discussed in Section II. 1, but with different initial condi- 
tions. During the expansion, the system gradually cools 
down and later decouples, when the temperature reaches 
a certain freeze-out value, Tj. In this model there is a 
simple relation between the temperature and the proper- 
time, i.e., r (X To • 

A very interesting alternative picture of the particle 
emission was proposed by Grassi, Hama and Kodama [sij: 
instead of emitting particles only when these crossed 
the freeze-out surface, they considered that the process 
could occur continuously during the whole history of the 
expanding volume, at different temperatures. In this 



model, due to the finite size and lifetime of the ther- 
malized matter, at any space-time point of the sys- 
tem, each particle would have a certain probability of 
not colliding anymore. So, the distribution function 
f{x,p) of the expanding system would have two com- 
ponents, one representing the portion already free and 
another corresponding to the part still interacting, i.e., 

f{x,p) = ffree{x,p) + fint{x,p). 

In the Continuous Emission (CE) model, the portion 
of free particles is considered to be a fraction of the to- 
tal distribution function, i.e., ffree{x,p) = Vf{x,p) 
fint{x,p) = (1 - V)f{x,p) or, equivalently, 

V 

ffree{x,p) = ^ _^ fint{x,p) ■ (58) 

The interacting part is assumed to be well represented 
by a thermal distribution function 

■f^'^'^'^^P^ ~ /t/.(a^,P) = (2^)3 exp [p.u{x)/T{x)] ± 1 ' 

(59) 

which poses a constraint on the applicability of the pic- 
ture, since by continuously emitting particles the system 
would be too dilute to be considered in thermal equilib- 
rium in its late stages of evolution. In Eq. (|59|l . is 
the fluid velocity at x'^ and T is its temperature in that 
point. 

The fraction of free particles V at each space-time 
point, x'^, was computed by using the Glauber for- 
mula, i.e., V = exp ^— n{x')avreidt'^ ,where tout = 

t + (-pcos0 -t- ^^Rl^- sin^ (/)) / (u sin 6*) . 

The model also considered that, initially, the energy 
density could be approximated by a constant (i.e., e — 
YqTq for all the points with p < Rt and zero for p > Rt)- 
Then, the probability V may be calculated analytically, 
resulting in V = {r/ToutY] a ~ 3^1^ Tg^Tocr-yrei, where 
Vrei ~ 1 • The factor V can be interpreted as the fraction 
of free particles with momentum p^ or, alternatively, as 
the probability that a particle with momentum es- 
capes from x^ without further collisions. 

Their early results for the spectra can be seen in 
Ref.[H3| and in the review by Frederique Grassi, in this 
volume, which is a good source of further details and 
discussions on the Continuous Emission model. 

Later, in collaboration with F. Grassi, Y. Hama, and 
O. Socolowski|53|. we developed the formulation for ap- 
plying this new freeze-out criterium into tttt interfer- 
omctry. Naturally, we would like to further explore if 
the above model would present striking differences when 
compared to the standard freeze-out picture (FO). One 
expectation would be that the space-time region from 
which the particles were emitted would be quite different 
in both scenarios. In particular, as we saw in Section II. 1 
and II. 2. 2, a non-instantaneous emission process strongly 
influences the behavior of the correlation function. Being 
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so, a sizable difference was expected when comparing the 
instant freeze-out hypothesis and the continuous emis- 
sion version, since in this last one, the emission process 
is expected to take much longer. 

For treating the identical particle correlation within 
the continuous emission picture, instead of using the Co- 
variant Current Ensemble formalism discussed in item 
II. 2, we adopted a different but equivalent form for ex- 
pressing the amplitudes in Eq. H26(l as in Kel. [5ll| . Then, 
the single-inclusive distribution, can be written as 



(60) 



neglecting the transverse expansion. For the instant FO 
case, we use Ea.(l5B)l and for obtaining C{q,K). For 
the CE case, we use Ea. (|5n)l and , writing the four- 
divergence and the integrals in cylindrical coordinates, 
using the symmetry of the problem, which leads to a 
simpler expression: 



G{q,K) 



(2^)^ Jo 



df] { pdp Tjr Mt cosh(y — if) 
Jo 



xe 



i[''':F(qo cos?) — g_L sinhjj)— pgr cos(0— </>,)] 



where 2?^ is the generalized divergence operator. In 
Ref.[E3|j it was shown that, in the limit of the usual 
freeze-out, Eq. (|50|) is reduced to the Cooper-Frye in- 
— Jrj.^dafj,p^f{x,p), over the freeze-out 



rdrpjrKT cos ip e 



'Mt cosh(y-?))/Tp^(x) 



,i[r((2o cosrj-gi, sinli rj) -pjcgr cos(0— 0,)] 



} 



(1 - V^ ) 



(63) 



tegral, 

hyper-surface T = Tf, being da^ the vector normal to where Mt = ^K'^ + K-P- = K^K^ - |g^(7,M = 

^{kj +k2); Y is the rapidity corresponding to K, (p is the 
azimuthal angle with respect to the direction of K, (j)q is 
the angle between the directions of q and K, and 



this surface. Or equivalently, Eq. (|60l ) in the instant 
freeze-out picture is reduced to Eq. with the cur- 

rents given by Eq. (|29(l . or even by Eq. (|3()|l . which is 
a simplified parametric form of describing thermal cur- 
rents used for obtaining analytical results in the Bjorken 
picture. We see more easily that it is indeed the limit 
if we replace the distribution function in Eq.|55J by its 
Boltzmann limit. We see more clearly that Eq. H28|l . or 
Eq. l|38(l in the Bjorken picture, are the natural limits 
of the proposed continuous emission spectrum, in case of 
instant freeze-out. 

Analogously, the two-particle complex amplitude is 
written as 



—p cos ^ -f -y/ J?|, — sin^ (p 



pT — —r (/ctZ-B) COS (;/)[a/ sinh^ + Pjr^''" — cosh ?;] ± 



K 



sinh^ ri + P. 



-2/a 



cosh ri] '■ 



G(fcl,fc2 



d^Xe''^^ {V^ W^ffree]}-- {2?^ Wiffree]}' 



(61) 

We saw above that the expression for the spectrum is 
reduced to the one in the Covariant Current Ensemble 
formalism in the limit of the instant freeze-out. Analo- 
gously, the above expression in this limit should yield to 
the result in Eq. H27|l . We see that this is indeed the 
case if we replace the individual momenta in E. H61|) 
by the average momentum of the pair, K^^ = 5(^1 + ^2 )• 
This replacement is even more natural, if we remember 
that the K^^ is the momentum appearing: in the Wigner 
formulation of interferomtrv fl^ Efll Is^l IIH . And, it 
was shown in Ref. lill that this also is reduced to the 
Covariant Current Ensemble for minimum packets and 
having the packet width equated to the pseudo-thermal 
temperature, as discussed in Section II. 2. 2. When this is 
assumed, also a substantial simplification is obtained in 
Ea. (|61|l . which could then be written as 



G{q,K) 



d^x e'""^- 



[K'^ffre 



(62) 



We compare next the results for these two very differ- 
ent scenarios by means of the two-pion correlation func- 
tions, assuming the Bjorken picture for the system, i.e.. 



The spectrum is obtained from the expression (|63|l . by 
replacing Rf" ^ k'^, M ^ m, and g'^ 0. In Eq. 
Tyr and pjt, whose expressions are written above, are the 
limiting values corresponding to the escape probability 
Vj^, which we fix to be Vj^ ~ 0.5, approximate value 
chosen for the sake of simplicity and for guaranteeing 
that the thermal assumption still holds in systems with 
finite size and finite lifetime. The rest of the emission for 
Vj^ > 0.5 is assumed to be instantaneous, as in Eq. 
and 

The complexity of the expressions for the Continuous 
Emission (CE) case suggested that we should look into 
special kinematical zones for investigating the differences 
between this scenario and the instant freeze-out picture. 
For details and discussions, see Ref.jH^- I summarize 
some of them. First, we observed that the correlation 
function C{qo,qL) plotted versus the outward momen- 
tum difference, qo, exhibited the well-know dependence 
on the average momentum of the pair, Kt^ discussed 
earlier in this section, in both cases. However, it was en- 
hanced in the CE scenario, as expected, since the emis- 
sion duration is longer in this case than in the standard 
freeze-out picture. We also observed a slight variation 
with Kt of the correlation function versus qs in the CE, 
differently from the instant freeze-out case, were it was 
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absent, since we were considering only the longitudinal 
expansion of the system and no transverse flow. This 
result showed the tendency of the curve versus qs to be- 
come slightly narrower for increasing Kt, an opposite 
tendency as compared to the curves plotted as functions 
of qo- 
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FIG. 18: The average correlation function {C)mv is shown 
versus ql, qo, andqs, averaging over, respectively, (qs , qo), 
(qh , qs), and (qh , qo)- Comparison is made of the (solid) 
curve for the CE scenario with three others corresponding to 
the usual freeze-out at Tfo = 140 MeV, but with different 
initial temperatures , wh ose values are shown in the plot. Plots 
extracted from Ref.\5^. 

We also studied more realistic situations, where the 
correlation function was averaged over kinematical zones. 
Using the azimuthal symmetry of the problem, we defined 
the transverse component Kt along the x-axis, such that 

K = [Kt, 0, Kl). We then averaged over different kine- 
matical regions, mimicking the experimental cuts, by in- 
tegrating over the components of q and K (except over 
the plotting component of q). For illustration, we con- 
sidered the kinematical range of the CERN /N A35 ,54J ex- 
periment on S+A collisions at 200 AGeV, as 



{C{qL))^l + 
ri80 r600 
-180 



/-ISO C dqs dq.C{K, q)\G{K, q)f 



/_i8o ho '^^^ Jo '^^s /o dloCiK, q)G{k^,h)G{k,, k^) 

With the previous equation we estimated the aver- 
age theoretical correlation functions versus ql, qo and 
qs, comparing the results for CE and for FO scenarios. 
First, we considered the case in which the initial temper- 
ature was the same in both cases (Tq = 200 MeV) and 
compared the CE prediction with two outcome curves 
corresponding to two FO temperatures, Tfo = 140 and 



Tfo ~ ITOMeV. We observed that the same initial tem- 
perature led to entirely different results in each case. The 
details are shown and discussed in Ref. js^]. In the sec- 
ond situation, we relaxed the initial temperature con- 
straint and studied the similarity of the curve in the CE 
scenario corresponding to a certain value of Tg, and com- 
pared to three different curves in the FO scenario, for 
which the freeze-out temperature was fixed Tfo = 140 
MeV. Each one of these curves in the standard FO case 
corresponded to a different initial temperature. The pur- 
pose here was to investigate, as usually done when trying 
to describe the experimental data, which initial tempera- 
tures in the FO scenario would lead to the curve closest to 
the one generated under the CE hypothesis. The results 
are shown in Fig. 18. Usually the shape of the correlation 
curve is very different in both cases, the one correspond- 
ing to the CE being highly non-Gaussian, mainly in the 
upper left plot of Fig. 18. In this particular one, we 
see that the CE correlation curve can be interpreted as 
showing the history of the hot expanding matter. For in- 
stance, the tail of (C) reflects essentially the early times, 
when the size of the system is small and the temperature 
is high, since the tail of the CE curve is closer to the 
FO one corresponding to the highest initial temperature. 
On the contrary, the peak region corresponds to the later 
times where the dimensions of the system are large and 
the temperature low (see the compatibility with the FO 
curve corresponding to the lowest initial temperature). 



II.6 FINITE SIZE EFFECTS 



II.6.1 PION SYSTEM 

In 1998, a pos-doctoral fellow from China, Qing-Hui 
Zhang, joined the group, spending one year with us. To- 
gether, we investigated the effects of finite source sizes 
(boundary effects) on tt^tt* interferometry. We derived 
a general formulation for spectrum and two-particle cor- 
relation, adopting a density matrix suitable for treating 
the charged pion cases in the non-relativistic limit [s^. 

A few hypotheses inspired that study. First, we con- 
sidered pions, the most abundant particles produced in 
relativistic heavy-ion collisions, to be quas i-bound in the 
system, with the surface tension |5^l5^l58l | acting as a re- 
flecting boundary. In this regard, the pion wave function 
could be assumed as vanishing outside this boundary. As 
usual, we also considered that these particles become free 
when their average separation is larger than their inter- 
action range and we assumed this transition to happen 
very rapidly, in such a way that the momentum distribu- 
tion of the pions would be governed by their momentum 
distribution just before they freeze out. We then studied 
the modifications on the observed pion momentum dis- 
tribution caused by the presence of this boundary. We 
also investigated its effects on the correlation function. 
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which is known to be sensitive to the geometrical size of 
the emission region as well as to the underlying dynam- 
ics. In this formalism, the single-pion distribution can 
then be written as 

Pi(p) = (vinp)^(p)> = EE'^A(p)v^v(p)(4«A') 

A A' 

= Y.Nxrx{p)Mp) , (64) 

A 

where the last equality follows from the fact that the ex- 
pectation value {a\a\') is related to the occupation prob- 
ability of the single-particle state A, N\, by {d\a\') = 
S\,\>N\ ; a\ {a\) is the annihilation (creation) opera- 
tor for destroying (creating) a pion in a quantum state 
characterized by a quantum number A. In Eq. H64|l 
^/;(t)(p) = ^r^V'/*'*(p)' where d'^ (d) is the pion cre- 
ation (annihilation) operator, and ipi is one of the eigen- 
functions belonging to a locahzed complete set, satisfying 
orthogonality and completeness relations. 

For a bosonic system in equilibrium at a temperature 
T and chemical potential /i, N\ it is represented by the 
Bose-Einstein distribution 

= TTTTT^ . (65) 

exp [^{Ex - fi)\-l 

The above formula coincides with the one employed in 
R,ef.[57l| for expressing the singie-pion distribution. 

The normalized expectation value of an observable A 
is given by 



Since we are considering the case of two indistin- 



(A) 



tr{pA} 
tr{p} 



exp 



T 



{H - pN) 



n 



pi 



(66) 

where p is the density matrix operator for the bosonic 
system, and pi = exp —^[Hi — /iiV;)j , with 

H = J2Hi; Hi= EiNi- N ^Y^^i ^ (^7) 

( I 

respectively, the Hamiltonian and number operators; p 
is the chemical potential, fixed to be zero in the results 
we present here. 

Similarly, the two-pion distribution function can be 
written as 

^2(Pl,P2) =E V'Ii(Pl)V'l2(P2)V'A3(Pl)'^A4(P2) 

Al,A2,A3,A4 

("Ai"A20A3aA4) 

E V'Ii(Pl)'/'l2(P2)V'A3(Pl)0A4(P2) 

Ai ,A2,A3,A4 



+ («Ai«A2^A3aA4)Ai=A2 = A3=A4 

A(Pi)Pi(P2) + |E^aV^a(Pi)Va(P2)|'. 



(68) 



guishable, identically charged pions, then {a\a}^axax) — 
2(d^dA)^. From the particular form proposed for the 
density matrix it follows that (dt^d^) = (axax) — 0, show- 
ing that it would not be suited for describing tt'^tt^ and 
7r"*'7r~ cases. For this purpose, the formalism proposed 
in Ref.jsil is more adequate. 

The two-particle correlation can be written as 



C2(Pl,P2) = 



= 1 



^2(Pl,P2) 
A(Pl)A(P2) 

IEa^a^I(Pi)^a(P2)P 



E,^a|^a(pi)PEa^^a|Va(p2)P ■ 

(69) 



In Ref. |55j we illustrated the formalism by means of 
two examples. The first considered that the produced 
pions were bounded inside a confining sphere with radius 
R. In the second, they were inside a cubic box with 
size L. Since the results were similar in both cases, I 
will briefly discuss here only the flrst one. We estimated 
Eq. (|64() for the confining sphere of radius R for studying 
the boundary effects on the spectrum. The results can 
be summarized by looking directly into the top left plot 
in Fig. 19, where we also show the curve corresponding 
to the limit of very large system (i? oo). We see 
that the finite size affects the spectrum by depleting the 
curve at small values of the pion momentum and, at the 
same time, rising and broadening the curve at large p 
(momentum conservation) . 

For studying the behavior of the correlation function 
we estimated Eq. in the case of the bounding sphere. 
Some results are shown in the top right plot of Fig. 19 
and they are as expected: the correlation shrinks (i.e., 
the probed source dimension increases) with increasing 
values of the radius. Nevertheless, the bottom curve, 
corresponding to C{q,K) versus q, shows an unexpected 
behavior for different values of the average momentum 
K. Contrary to what is observed in expanding system, 
the correlation function becomes narrower (probed re- 
gion enlarges) with increasing K. However, the exam- 
ple shown for illustration does not take into account the 
expansion of the system. The variation with K merely 
reflects the strong sensitivity of the results to the dynam- 
ical matrix adopted in the formulation. It is related to 
the weight factor Nx, in Eq. H69() : pions with larger mo- 
menta come from higher quantum states A, which cor- 
respond to a smaller spread in coordinate space. But, 
due to the Bose-Einstein form of the weight factor, large 
quantum states give a small contribution to the source 
distribution, causing the effective source radius to appear 
larger. To emphasize this we include in the bottom plot 
the narrowest curve, corresponding to flxing Nx = 1 (or 
equivalently, by considering T 1 in the Bose-Einstein 
distribution, which makes C{q,K) insensitive to K, due 
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Eq. (70) 




\ o numerical results 







q(GeV/c) 

FIG. 19: The top plot shows the spectrum in the example 
of pions confined in the sphere of radius R. The top right 
plot shows CTTTviq) versus q, for three different values of the 
source radius. The plot in the bottom shows that the correla- 
tion function shrinks for increasing average momentum, K , 
in opposition to the behavior seen in Fig. 7 and for expanding 
systems, in general. Plots extracted from Ref.JSSf. 



to the large temperature). This limit allows for an ana- 
lytical expression, 



C2{q) = l + - 



Rcos{qR) 



sm{qR) 



q 



(70) 



which is shown in the lowest curve of the bottom plot of 
Fig. 19. Underneath this curve, shown by little circles 
in the same plot, is the numerically generated curve for 
iV^ = 1, confirming the correcteness of our result. 



II.6.2 Q-BOSON SYSTEM 

More recently, Qing-Hui Zhang and I extended the 
above formalism for treating the interferometry of two 
Q-bosons. The concept of quons was suggestedje^ as an 
artifact for reducing the complexity of interacting sys- 
tems, at the expense of deforming their commutation re- 
lations by means of a deformation parameter, Q. Then, 
this could be seen as an effective parameter, encapsulat- 
ing the essential features of complex dynamics. 

We derived a formalism suitable for describing spec- 
tra and two-particle correlation function of charged Q- 
bosons, which we considered as bounded in a finite vol- 
ume. We adopted the Q-boson type suggested by An- 
chishkin et al. 61], according to which the Q-bosons 
are defined by the following algebra of creation {a'') 
and annihilation operators (a): aiaj, — Q^'^'' aj,di = 



it 



',] = 0, [Ni,ai,] = -Sij,ai [Ni, 



- 0; Ni is the number operator, Ni — 
'{aiY. In the limit Q = 1, the regular 



diMaJ, [Ni,Ni,] 

sr°° (1-9)° i^^t^ 
i^s=i (i-Q=) > 

bosonic commutation relations are recovered. The defor- 
mation parameter Q is a C-number, here assumed to be 
within the interval [0, 1]. 

The single-inclusive, -Pi(p), and the two-particle dis- 
tributions are derived in a similar way as in the case 
of regular pions discussed above. In the Q-boson case, 
Ea. (|64() continues to hold, but the weight factor, N\, re- 
lated to the occupation probability of a single-particle 
state no longer is as written in Ea. H65|l . but is changed 
into a modified Bose-Einstein distribution. 



1 



exp 



m)] -Q 



(71) 



The two-particle correlation function as also modified, 
and is written as 



C2(Pl,P2) 



-P2(Pl,P2) 



Pi(pi)Pi(p2) 

i + |^7v,i^,(pi)PE^'i^'(p2)i'j 



Ll 



exp(i(i;,-/i))-Q2 



(72) 



We apply this formulation by means of similar models 
as in Ref.|55j, for the confining sphere of radius R. We 
studied the effects of different values of Q on the correla- 
tion function, under different values of the pair average 
momentum, Kt. In Fig. 20 we show the boundary ef- 
fects on the Q-boson spectrum and on the two-Q-boson 
correlation function. We also included in that plot the 
behavior of the of the two-Q boson intercept parameter, 
A, defined by A = C(q = 0, K) — 1. The top two plots in 
Fig. 20 reproduce the basic characteristics seen in Fig. 
19, in particular, the correlation function narrows as the 
average momentum, K, grows. There is, however, a new 
result: the maximum of the correlation function drops as 
the deformation parameter, Q, increases. This can also 
be seen by looking into the bottom plot of the intercept 
A, as a function of Q. 

We also derived a generalized Wigner function for the 
Q-boson interferometry, which would reduce to the reg- 
ular one in the limit of Q — s- 1. For that, we define the 
Wigner function associated to the state / as 



ffi(x,K) 



d^Ax 
^2^ 



iK.Ax 



Ax, , , Ax, 

^r(x+— )v^((x-— ) 



(73) 

We proceeded analogously to define the equivalent 
function for the integration in y and Ay, remembering 
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- K = 0.3 GeV/c, R = 3 fm 
K = 0.5 GeV/c, R = 3 fm 
a K = 0.3 GeV/c, R = 6 fm 
^--o K = 0.5 GeV/c, R.6fm 



FIG. 20; The results shown here were obtained by consider- 
ing the Q-bosons confines inside a three-dimensional sphere 
of radius R. In the three plots above, the temperature was 
fixed to r = 0.12 GeV and the chemical potential, /i = 0. The 
top left plot corresponds to the normalized spectrum (in ar- 
bitrary units) vs. momentum |p| (in GeV/c). The solid lines 
correspond to the case i? — > oo, the dotted ones to i? = 6 
fm, and the dashed lines correspond to the case, R = 3 fm/c. 
The bare lines refer to Q = 1.0 and the ones with symbols, 
to Q = 0.5 ((3 = 0). The top right plot shows the two-pion 
correlation vs. momentum difference |q| (in GeV/c), for a 
i? = 3 fm sphere. The solid lines correspond to mean mo- 
mentum K — 0.3 GeV/c and the dashed ones, to the case 
if = 0.5 GeV/c. The bare lines refer to the case of Q = 1.0 
and the ones with symbol, to Q = 0.5 (Q = 0). The bottom 
plot shows the intercept parameter. A, vs. the deformation Q. 
The solid lines correspond to the mean momentum if = 0.3 
GeV/c and dashed ones to the case if = 0.5 GeV/c. The bare 
lines refer to the case i? = 3 fm and the ones with symbols, 
to i? = 6 fm. Plots extracted from Ref . [62^] . 



that gi{x,Ji) = g,*(x, K). Then, dcnotmg by 

g{^,K)=Y.^i9i, (74) 
I 

we finaUy defined the generalized Wigner function of the 
problem as 

5(x,K;y,K) =5(x,K)g(y,K) - (1 - Q) x 



exp(i.(i;i-M))-Q2 



5/(x,K).9,(y,K) .(75) 



We see that, for Q = 1, the above expression is reduced 
to the usual result of the original Wigner function, i.e., 



g(x,K;y,K) =5(x,K)g(y,K) 



(76) 



On the other hand, for Q = Q, Eq. H75|l is identically zero 
for single modes only, as it would be expect in the limit 



of Boltzmann statistics. Nevertheless, in the multi-mode 
case, there seems to be some sort of residual correlation 
among Q-bosons even in the classical limit. 

By means of this Wigner function, the two-Q-boson 
correlation function can be rewritten as 



C2(Pl,P2) = 1 



-4q.(x- 



-y)g(x,K;y,K)dxdy 



/ Pi)dxJ g{y, P2)dy 



(77) 



Particularly interesting is the approach in Ref. |63| . 
where it was shown that the composite nature of the 
particles (pseudo-scalar mesons) under study could re- 
sult into deformed structures linked to the deformation 
parameter Q. In that reference this parameter is then 
interpreted as a measure of effects coming from the in- 
ternal degrees of freedom of composite particles (mesons, 
in our case), being the value of Q dependent on the de- 
gree of overlap of the extended structure of the particles 
in the medium. Being so, the Q-parameter could be re- 
lated to the power of probing lenses, for mimicking the 
effects of internal constituents of the bosons. In this case, 
and for high enough magnification, the bosonic behavior 
of the Q-bosons could be blurred by the fermionic effect 
of their internal constituents, which would result in de- 
creasing the value of Q. Our results on the two-Q-boson 
interferometry are compatible with this interpretation, 
as explained in detail in Ref.(22i- 



II.7 NON-EXTENSIVE STATISTICS AND HBT 

Sergio M. Antunes[64l|. working under my supervision 
and in collaboration with G. Krein during his Master 
Degree, studied the effects of Tsallisj^^ non-extensive 
statistics on the two-particle correlations and spectra. In 
this work, a very simple starting hypothesis was made: 
that under certain circumstances, the Boltzmann limit to 
the pion distribution could be replaced by an equivalent 
approximate expression derived within the non-extensive 
statistics. This concept of non-extensivity could be sum- 
marized very briefly by the relation: {A + B) = 
Sg^^iA) + Sg^^iB) + (1 - qTjSg^^iA)Sg^jB), i.c., the 
entropy of a system formed by two independent sub- 
systems A and B, no longer is the sum of the entropy 
of the two subsystems (note that by independent it is 
meant that the probability of the composite system fac- 
torizes as pa+b = PA + Pb)- The parameter qxs is a 
measure of the degree of non-extensivity of the system 
and the Boltzmann-Gibbs statistics is recovered in the 
limit qx^ — > 1. From the definition of the generalized 
entropy in the Tsallis.65.j formulation, it is possible to 
obtain approximate analytical expressions, for instance, 
for the mean occupation number |66l| 



exp{h(£;(p,r)-/i]/r}^ 

{l + (gT. -l)[(i?(p,r)-/.]/^}-l/^^--- 



(78) 



22 



which is reduced to the Boltzmann distribution for qx^-^l, 
where T is the temperature and fj,, the chemical potential. 
The form given in Ea. l|78() is, however, valid for qxs close 
to unity. Later, G. Wilk et al.[63 showed that the above 
distribution can be written in the form 



Gn 



Co 



1-(1-9tJ 



A 



(79) 



Considering Gq^^ as a probability distribution (Levy dis- 
tribution) in the variable x, with x e (0,cxd), the param- 
eter qTs must be limited to the interval 1 < qTs < 2. 
However, if the mean value of x is required to be finite 
(< X >= A/(3 — 2qTs) < oo) for x e (0,oo), then q^s 
has to belong to the interval 1 < qTs < 1-5, which better 
justifies the interval of applicability of Eq. H78|) . 

For investigating the above hypothesis, we adopted a 
model with radial flow 68], leading to a non-decoupled 
phase-space freeze-out distribution 

''^"'^^^[l + (.T.-t^^}^,^i-M)]^/(--'' ^''^ 

where E'{p,f) — {Ep — r.p)(l — u^)~^ and /i = 0. 

With the above decoupling distribution the correlation 
function is written as 

C2{K,q) = l + exp[-^{Ek, - Eu^fr^] x 



pTT sin 9jQ{\q \ |:r jsin 9 sin 6) exp{2|(7| |:r|cos 9 cos 9^d6 



[l + (qT,~l){i:E 



T 



-J/ cose)] 



l/(9Ts -1) 



{l^' d([l + {qn - l)(^^i?/2 - ycos0)]i/(i-.) 



where Ei ^ 



(f )2 m2 , 7 = (1 - «2)-l/2 j^(^) 

is the Bessel function of order 0. In Eq. 1)81(1 . 6 is the 
angle between K and x, 9 is the angle between K and 
g, in such a way that the angle a between q and x, can 
be determined by cos a = cosOcosO -\- sinOsinO cos{4> — 
(j)). But K and q define a plane, so that we can choose 
the directions of these vectors in such a way that 4> = 
0. Then, choosing the direction of K along the z-axis, 
if we integrate Eq. (jHU for ^ = and ^ = 90°, this 
corresponds, respectively, to q\\ K (i.e, qo) and q Jl K 
(i.e., qs). 

The correlation function versus the momentum differ- 
ence, q — ki — k2, showed a strong dependence on the 
combined variables K — ki + k2 and T/(7v), where v is 
the flow velocity (with v ~ 0.55 as typical flow velocity 
at CERN/SPS, and v ^ 0.37 as typical flow velocity at 
LBNL/Bevalac), 7 = (1 — f^)"^, and T is the tempera- 
ture. Although not explicitly illustrated here, the depen- 
dence on the average momentum of the pair showed that 
the correlation curves became narrower for increasing K , 
similarly to what was observed in the results of Sect. II.6, 
even though the present model considers an expanding 




FIG. 21: The two upper plots show results for the two-pion 
correlation as function of qo- The top plot has fixed values 
of l^j = 0.8 GeV/c and the factor T/i^yv) = 0.152 GeV, with 
T = 0.100 MeV and v ^< v >— 0.55, and curves correspond- 
ing to ~ 1.0 (dashed), qTg ~ 1.015 (long-dashed) and 
QTg =1.10 (continuous). The plot in the middle results from 
fixing 1^1 = 0.6 GeV/c, T/{jv) = 0.109 GeV, T = 0.044 MeV 



and V «< V >= 0.37, with curves corresponding to qr^ 



1.0 



(long-dashed) and qxs = 1.015 (continuous). The plot on 
the bottom show the spectrum corresponding to the upper 
case, with corresponding curves for gyg = 1.0 (continuous) 
and qTs ~ 1.015 (dashed). 



system. Moreover, also similarly to the results on the 
Q-boson interferometry, the correlation shrinks, i.e., the 
probed effective region grows, for increasing qx^ , suggest- 
ing that long-range correlation could be present, in asso- 
ciation with Tsallis statistics. The study also showed that 
a very small deviation from the Boltzmann statistics, cor- 
responding to qT^ — 1.015, could lead to clear differences 
in the correlation function under some particular choices 
of the combination T/{'jv), as shown in the top right 
plot of Fig. 21. We see that, due the strong dependence 
on the combination T/^jv), the search for such a devia- 
tion from the Boltzmann statistics, as suggested by the 
Tsallis extensive statistics, would be favored at lower en- 
ergies. Nevertheless, it was shown that the experimental 
data on event-by-event pion transverse-momentum fluc- 
tuations, that could not be explained by a model based 
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on standard extensive statistics, was compatible with a 
small deviation, i.e., for a value of the non-extensive pa- 
rameter qx^ = 1.015, which inspired our analysis. Also, 
NA35 Collaboration data on pion transverse momentum 
distribution from S + S collisions at SPS showed a better 
agreement with fit based on non-extensive statistics, with 
qx^ = 1.038|73, mainly for the tail of the distribution, 
which is of power-law type. 

The above approach to the problem, however, is very 
simple, since the connection of qr^ to the space-time 
variables was considered only through the non-decoupled 
phase-space present in the adopted modeljH^- Never- 
theless, the possibility is not excluded that a more gen- 
eral form for the Wigner function, as the one derived 
in Ea. (f75|l of Section II. 6. 2, would be more appropriate. 
Maybe long-range correlations, as suggested by Tsallis 
statistics, would not allow for reducing the Wigner func- 
tion to its conventional form, as in Ea. (|76|) . 



II.8 PHASE-SPACE DENSITY 

Victor Vizcarra-RuizlXll, another student working 
with me, investigated a suggestion made by George 
Bertsch|72, according to which the average phase-space 
density could be estimated by means of the two-pion cor- 
relation function and the single-particle spectrum. In his 
Master dissertation, a generalization of Bertsch's sugges- 
tion is proposed, by applying the wave-packet formalism 
proposed in Ref.|4l|. 

Bertsch's proposition could be briefly summarized as 
follows. He considered that, in ultra-relativistic heavy- 
ion collisions, the local phase-space density in the final 
state is frozen and gives a measure of the dynamics in 
the priory interacting region. He starts by converting 
the source function to an equivalent one at a common 
instant, to, times the phase-space density at that time, 
i.e.. 



gix,K)^S{t-to)f{r,K)/{27r)' 



(81) 



Using the Wigner formulation for the spectrum = 
/ cPxg{x, K) and assuming the approximate validity of 
the relation dPN/d^ki « (fN/d^k2 « cfN/tfiK, inde- 
pendently of the momentum difference q, he wrote the 
average phase-space density as 



<f>K- 



Jd^rfir^K) _ d^N , 



Jd3rfif,K) 



d^k 



A[C(g,X)-l]. 



(82) 

For instance, in case the experimental single-inclusive 
distribution could be well reproduced by the expression 



d^N _ 1 dNy e-^-rlT 



(83) 



and the correlation function by Ea. (|17|l . the phase-space 
density would be written as 



A0F dNy e-^^/'^ 
2EkT^ dy RoRsRl 



(84) 



since the right-hand-side of Eq. H82|l results in ^^^^^^ . 

The generalization proposed in |7l| was derived us- 
ing the powerful formulation based on the Wigner 
formalism |4l| and summarized in Section II. 2. 2, applying 
to the case of two-boson interferometry. In this first ap- 
proach, we tried to keep our derivation as close as possible 
to the one proposed by Bertsch, so that the differences 
could be clearly seen. With this in mind, we maintained 
expression (I81|l for the non-relativistic source function as 
suggested in |7a |. 

The single-particle distribution (spectrum) is then 
written as 



d^N 



d^h 



1 



(27r)3 

where we have defined 
f{f,p,K;Ap) = 



A (27r)3(27rAp)t 

d^rd^p f{f,p,K;Ap), 



d\d^pf{r,p) e-(P-^»)V2Ap^ 



(85) 



1 



(27rAp)t 

The average phase-space density is then defined as 



(86) 



< f{K; Ap, Ax) >= 

J d^rd^p f{r, p, K; Ap, Ax) f(r, p, K; Ap) 
J d^rd^p f{f, p, K; Ap) 



(87) 



Then, the reformulation of the average phase space 
density in terms of the generalized Wigner formalism was 
derived as 



<f{K;Ap,Ax) 



d^N 



d^K 



d\[CA{q,K)-l] 



With the definition in Ea. l|HHjl . we see from Eq. H88(l 
that the form of expression for the average phase density 
< / > does not change significantly with the introduction 
of the wave-packets coming from the generalization of the 
Wigner formulation for interferometry,. What actually 
changes is the way to estimate the spectrum and the 
two-pion correlation function, where the wave packets 
are implicitly present. 

For illustrating an application of Ea. (|88|l . we can look 
into the simplest example, the Gaussian emission func- 
tion, similarly to the one proposed by Bertsch, but with- 
out splitting the transverse variables, i.e., 



C{ki,k2) = 1+ e 



— qrpRrp —q^R L 



(89) 
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FIG. 22: The top plot shows the pion spectral distribution. 
The dashed line represents the transverse distribution in Eq. 
H8i-i|l . the continuous one, the Kolehmainen-Gyulassy type, 
and the long-dashed curve, the Gaussian spectral distribution. 
The plot in the middle shows the correlation function in terms 
of ql, for all the parameterizations used: Gaussian without 
and with the inclusion of the wave-packets (two bottom ones) , 
as well as the one on top corresponding to minimal packets (or 
Kolehmainen-Gyulassy). The plot on the bottom shows the 
estimate for the average phase-space (divided by the integral 
in rapidity), for all the models: the highest curve corresponds 
to the simple Gaussian parametrization for the source with 
exponential spectrum, the lowest one, to the Gaussian source 
with Gaussian spectrum and minimal wave-packets, and the 
middle curve, to the Kolehmainen-Gyulassy type. 



whereas the correspondent exp ression in the generaUzed 
Wigner formulation of Ref . |41| is 



C(fcl,fc2 



1 



(90) 



from where we see that the intercept parameter, A, is im- 
pHcitly being fixed to unity. In Eq. , we have R\ ^ ~ 
RfiK) + l/(4m^T/). Thus, we see that the results for 
the radial parameters in the standard Wigner formula- 
tion would correspond to the hmit Ax ~ 1/Ap — 0. 
The parameters used in this example are: Rt = 3.56 fm, 
Rl = 2.52 fm, as suggested by experimental results from 



the BNL/AGS|4g, Tf = 140 MeV, and minimum packets 
were considered, i.e., AxAp = 1/2, being Ap w yJmT^T^. 
Relative to these values of radii, the factor 1/Ap is very 
small (0.5 fm, or « 4%). Consequently, the corrections to 
the correlation function and to < / > are expected to be 
small, which is indeed what is observed. Nevertheless, 
we should notice that the time dependence was supposed 
to be a delta function, similar to Eq. (|^ . to keep the 
analysis simple and closer to Bertsch's proposition, and 
for emphasizing the influence of the wave-packets. How- 
ever, we know that the duration of the emission has an 
enormous influence into the correlation function and the 
expectations are that more pronounced differences would 
be seen if a broader time interval was introduced, for in- 
stance, by Gaussian time spread instead of the delta one 
used. 

The proposed extension of the average phase-space es- 
timate was applied to some examples. The first was the 
static Gaussian case, with and without the inclusion of 
wave packets. Another one was the case corresponding 
to the Bjorken picture, IOC, which, for minimal pack- 
ets in the generalized Wigner formalism, is reduced to 
the Kolehmainen-Gyulassy solution, where the spectrum 
is close to an exponential. In the generalized Wigner 
formalism, however, it is possible to choose other possi- 
bilities for parameterizing the spectrum and, for a test, 
a Gaussian version of the momenta distribution is tried, 
although it is well known that the experimental results 
favor an exponential type of parametrization. 

In the estimate of the average phase-space density, a 
strong sensitivity to the adopted spectral shape was ob- 
served. The study compared cases of static Gaussian 
source distribution, having exponential and Gaussian 
corresponding to minimum packets, with the Bjorken 
type from the Kolehmainen-Gyulassy approach. This is 
shown in Fig. 22. 



II. 9 THE RHIC PUZZLE 

The new millennium started together with the first 
run of the Relativistic Heavy Ion Collider, at BNL. The 
main preliminary surprising result released was, coinci- 
dentally, in HBT. Since The conference Quark Matter 
2001, in January of that year, it has been known as the 
RHIC Puzzle, first seen by the STAR collaboration!?^ 
and later confirmed by PHENIX Collaborationlz^- The 
result was related to the ratio of two of the radii param- 
eters, Rout/Rsid (Rout = Ro is the transverse radius in 
the direction of the transverse average momentum, Kt, 
whereas Rsid = Rs is its component orthogonal to Kt). 
The predictions based on microscopic models, such as 
UrQMD, suggested that the ratio would either grow in- 
definitely (considering hadronic re-scattering only) or in- 
crease up to Kt > 0.15 GeV/c and then decrease, appar- 
ently saturating towards unity at large Kt (when com- 



25 



billing hydrodynamics and cascading). In any case, the 
ratio was above unity. The data, however, indicated an 
opposite trend, i.e., the data points for that ratio ap- 
peared to be steadily decreasing, reaching values around 
0.8 at Kt « 0.4 GeV/c. Hydrodynamical models, such as 
proposed by U. Heinz and P. Kolbj75||, failed in general to 
describe the experimental results for the individual radii, 
approaching data only if they considered that the sys- 
tem would immediately freeze-out at the hadronization 
point. Even in that case, however, they did not succeed 
in describing the ratio. 



II.9.1 HOT TAMALE MODEL 

Challenged as well by this puzzle, Larry MacLerran 
and myself tried and built a model assuming that the 
QGP, initially formed at RHIC, constituted an opaque 
source. The opacity was implemented by considering a 
model where pions are emitted from the surface of the 
system, at fixed radius, all along its lifetime, from its for- 
mation at To, up to the freeze-out. We neglected trans- 
verse flow for simplicity, considering the ideal Bjorken 
picture of 1+1 longitudinal expansion. The phase tran- 
sition started at Tc (at temperature Tc = 175 MeV, to be 
consistent with lattice Monte-Carlo data), ending at r^,. 
The system, pions only, further expanded until reaching 
Tf (at Tf = 150 MeV, to be consistent with the energy 
per particle at RHIC and to fit the pr distribution of 
pions), when it broke up, in a volumetric emission. 

Many of the features of the model we propose are em- 
bodied in the hydrodynamic computations of Heinz and 
Kolb[l3, the essential difference being the treatment of 
surface emission. The energy was emitted already from 
the partonic phase. This quark and gluon matter was as- 
sumed to be directly converted into a flux of pions with 
the same energy and a blackbody distribution at the tem- 
perature of emission. This energy conservation condition 
allowed us to directly take a fiux of gluons and quarks and 
convert it into a spectrum of purely pions. No complex 
mechanism for the QGP hadronization was considered in 
detail, although hadronization must take place. In other 
words, in first approximation, we considered the evap- 
oration of "gluons" and "quarks" (as hadronized pions) 
from the external surface of the system in the same way 
as emission of pions, except for the different number of 
degrees of freedom. 

We estimated the emitted energy and the total entropy 
at each stage. In the initial phase, lasting from tq to 
Tc, the emitted energy as a function of time was esti- 
mated considering the emission by an expanding cylinder 
of transverse radius Rt and length h, in the time interval 
between r and t + dr, as 



^hkjT'^ 2t: Rt hdr - \ aT'^ tt R^dh , (91) 



where the first term comes from the black-body type of 
energy radiated from the surface of the cylinder, and the 
second term results from the mechanical work due to its 
expansion. The k factor was introduced to take into ac- 
count that the system has some opacity to surface emis- 
sion. The constant a is proportional to the number of 
degrees of freedom in the system. Integrating this equa- 
tion we get for the energy density (i.e., e = E/V) 



(92) 



Em = eo( — )3 e °' . 

r 

From Ea. (|92|l we see that an extra factor, 
obtained in additional to that coming from the Bjorken 
picture. Consequently, the temperature in this model 
changes more rapidly, according to 

T(r) = ro(-)4e"™?^"""°^ (To<T<T,) . (93) 

T 

The initial temperature was estimated by equating 
the initial entropy to the number of produced pi- 
ons (pion yield) at the end, i.e., = FTV = 
(5s + 9q) X (|)m^o 7ri?|, To, where M ~ 1000 is the av- 
erage produced pion multiplicity per unit of rapidity at 
RHIC, and F = S^/N^ « 3.6. Then, we get To = 411 
MeV. The degeneracy factors, g, are given by the gluon 
degrees of freedom, gg = 2(spin) x 8(color), and the 
quark/anti-quark degrees of freedom, gg = |[2(spzn) x 
2(9 + 9) X 3{color) X Nf {flavor)], which add up to 
9qgp — 9g + gq- In the case of pions, the degeneracy 
factor is = 3. 

The initial time is estimated by means of the uncer- 
tainty principle, i.e., EqTq ^ 1 [he = 1) and Eq « 3ro. 

TABLE 4: Parameters of the Tamate Model 



K 


To 

{¥) 


Tc 
(^) 


Th 
(^) 


Tf 
(^) 


S/Aftot 

(to < T < Tf) 


V/^tot 

(at Tf) 


1 


0.160 


1.54 


5.73 


6.97 


0.844 


0.156 


0.5 


0.160 


1.75 


8.37 


10.5 


0.758 


0.242 



During the mixed phase (i.e., Tc < t < Th), the tem- 
perature remains constant and the total entropy is con- 
served. If the fraction of the fluid in the QGP phase is 
/, the fraction in the pionic phase is (1 — /), so that the 
entropy density is written as s = fsqcp + (1 — f)sh, 
where 



/ 



T e 



Th ■ 



Tr e 



-Ath- 



(94) 



Finally, the instants corresponding to the initial of the 
first-order phase transition, t^, its end, Th, and the 
breakup of the system, Tf, were estimated as 



To e- 



(95) 
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Table 4 illustrates the time variables for two different as- 
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FIG. 23: Part (a) shows the prediction for the spectrum based 
on the Tamale model. The points are from the minimum 
bias data from PHENIX Collaboration. The emissivity k is 
indicated in each curve. Part (b) shows results for the ratio 
Rout/Rsid for our model and for two values of the emissivity. 
The ratio corresponding to k = 1 agrees very well with the 
PHENIX data points, within the experimental error bars but 
the one with 50% emissivity is completely excluded by data. 
The plots were extracted from Ref . |7al . 



sumptions on the emissivity, n. I also write the estimated 
fraction of the particles emitted from the surface during 
the period tq < t < r/, S /Ntot, relative to the total num- 
ber of produced particles, Ntot, as well as the remnant 
portion at freeze-out, "V /Mtot, then emitted from the en- 
tire volume. With the above parameters, we estimated 
the spectra and compared with experimental data from 
RHIC/Phenix CoUaborationlll, restricted to the central 
rapidity region. This is shown in Fig. 23. 

With this model we managed to roughly explain the 
trend of data for Rout but failed for Rsid, since it was 
independent on Kt- This, however, was expected, since 
the flow was neglected. The ratio, however, as can be 
seen from Fig. 23, came up just right, through the mid- 



dle of the data points^]! Nevertheless, only much after- 
wards we did realize our emissivity parameter, k, which 
controlled the surface emission, differed from the conven- 
tional one associated to a black-body, fn fact, only by 
then we realized that the system would have to irradi- 
ate about four times more than an usual black-body for 
matching the data points. We then concluded that it was 
another strong evidence that flow should be included in 
the estimates, which is still under investigation. 



II.9.2 IC, SPheRio, CEM & HBT 

In Ref. [ill J Otavio, Frederique, Yogiro and Takeshi an- 
alyze the effects of fluctuations in the initial conditions 
(IC), and of the continuous emission (CE) model, in two- 
pion interferometry within the hydrodynamical descrip- 
tion of the SPheRio codejl^, comparing the results with 
the RHIC/STAR data. 

The IC are, for simplicity, usually chosen as highly 
symmetric and smooth distributions of velocities and 
thermodynamical quantities. These IC correspond to 
mean distributions of hydrodynamical variables, aver- 
aged over several events. Nevertheless, for the typical fi- 
nite systems formed in high energy collisions, large fluctu- 
ations are expected, varying from event to event. More- 
over, the IC on the event-by-event basis often show small 
high-density spotsjT^ in the energy distribution. Being 
so, it would be expected that such spots would manifest 
themselves in the particle emission, contributing to de- 
crease the HBT radii. In order to produce event-by-event 
fluctuatingIC, in Ref. [ill they used the NeXus event 
generator [79j: once the incident nuclei and incident en- 
ergy are given, it produces the energy-momentum tensor 
distribution at the proper time ^t"^ — z^, in an event- 
by-event basis. This, together with the baryon-number 
density distributions constitute the fluctuating IC. No 
strangeness was introduced in the calculation. 

On the other hand, the CE model we already discussed 
in Sec. II. 5, is a alternative picture for describing the fi- 
nal state of the high energy collisions, i.e., the particle 
emission. According to the CE model, this emission oc- 
curs not only from the sharply defined freeze-out surface 
but also happens continuously, from the whole expanding 
volume at different temperatures and times. The main 
ingredients of the model are discussed in more detail by 
Frederique Grassi in her review, in this volume. The ba- 
sic formulation for CE interferometry was discussed in 
Sec. II.5. 

In between the assumptions in the IC and the final 
stages described by the CE picture, there is the system 
evolution. This is considered in Ref.[73 by means of hy- 
drodynamics, developed in the code SPheRio, which was 
based on the Smoothed Particle Hydrodynamics (SPH), 
first used in astrophysics and more recently adapted to 
nuclear collisions [T^. SPH uses discrete Lagrangian co- 
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ordinates attached to small volumes ("particles") with 
some conserved quantities, which they take as the en- 
tropy and baryon number. For details about SPheRio 
and hydrodynamics, see the review by Yogiro Hama and 
Takeshi Kodama, in this volume. 

Some of the results shown in Ref . (77| are summarized 
in Fig. 24. To test the effect of fluctuations in the IC, 
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FIG. 24: The plot on the left shows correlation functions 
assuming fluctuating IC and average IC, with sudden freeze- 
out. The qo,s,L that are not displayed in the horizontal axes 
were integrated in the interval < qo,s,L < 35 MeV. The 
plot on the right shows the radii Ro,Rs,Rl together with 
the ratio Ro/Rs, as a function of ttit. The complement 1 
stands for average IC and 2, for fluctuating IC. The plots 
were extracted from Ref. [t^ . 

a sudden freeze-out is initially assumed, at T^.o. — 128 
MeV. In the left part of Fig. 24 the correlation function 
C2, averaged over 15 fluctuating events, is compared with 
the corresponding C2 with average IC (i.e., no fluctua- 
tions). The IC fluctuations are reflected in the fluctu- 
ations seen in the correlation functions. Averaging the 
correlations over the fluctuating range, we see that the 
resulting C2 is broader than the correspondent ones com- 
puted with average IC. So, the effect of the IC fluctua- 
tions is to reduce the effective radii extracted from the 
correlation functions. Besides, the shapes of the result- 
ing correlations with fluctuating IC and with smoothed 
average IC, are different. On the right-hand-side of Fig. 
24, it is shown the tot dependence of the HBT radii. 
These were estimated by fitting Gaussians to C2. Re- 
sults based on the CE model, as well as RHIC data from 
STAR Collaboration 'tsJ are also included in that plot. 

From Fig. 24, we see that joint effects of the smooth IC 
with sudden freeze-out make the uit dependence of the 
outwards radii, Rq, flat or slightly increasing. The fluc- 



tuating IC push the effective radii towards smaller val- 
ues and flat niT dependence. Adding the CE hypothesis 
to the analysis, we also see from Fig. 24 that all effec- 
tive radii, Rq, Rs, and Rl decrease with increasing rriT. 
This is expected, since in the CE picture, high-fc-r parti- 
cles have a bigger chance of being emitted at early times, 
when the system is hot, mostly from its surface. On the 
contrary, low-kT particles are emitted later, when the 
system has already cooled down and has larger size due 
to the expansion. The role of the fluctuating IC in the 
CE case is mostly of reducing the values of the effective 
radii, without significantly changing its uit dependence. 
These two effects together improve the description of the 
data points, as seen in Fig. 24. Nevertheless, the ratio 
Ro I Rs is still above unity and does not describe the data 
points yet. Another ingredient may still be lacking. 



11.10 BACK-TO-BACK CORRELATIONS 

The Back-to-Back Correlations (BBC) have a much 
more recent history. Under the historical perspective we 
are pursuing in this review, it is worthy to briefly revise 
its origin. In 1991, Andreev, Pliimber and Weinerjsojl 
wrote a paper in which they pointed out the surprise ex- 
istence of a new quantum statistical correlation, among 
particle- antiparticle pairs. Thus, tt+tt" correlations of 
this type would be similar to tt'^tt*' (since 7r° is its 
own antiparticle) but entirely different from the Bose- 
Einstein correlations we discussed in the previous sec- 
tions. They were related to the expectation value of 
the annihilation (creator) operator, more specifically, 
< a('l')(fci)a(''^^(fc2) >7^ 0, analogous to what is observed 
in two-particle squeezing in optics, where the averages are 
estimated using a density matrix that contains squeezed 
states. Although not all they discussed was correct, they 
pointed out that C{tt+tt^) > 1 and C(7r''7r'') > 1, re- 
flecting particle- antiparticle quantum statistical effects. 
Later, SinvukoV|i59t|. discussed a similar effect for tt+tt" 
and 7r°7r" pairs, claiming that the effects would be due 
to inhomogeneities in the system, as opposite homogene- 
ity regions in HBT, which comes from a hydrodynamical 
description of the system evolution. He used Wick's theo- 
rem for expanding the two-particle inclusive distribution 
in terms of bilinear forms. 

We can understand the origin of this effect in a simpler 
way in terms of creation and annihilation operators, in 
the tt'^ case. For instance, the single- inclusive distribution 
is written as 



A^i(k. 



dk,. 



and, after the decomposition that follows from Wick's 
theorem, the two-particle distribution is written as 



Af2(ki,k2) = t^ki^ka 
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WkiWki{(4iaki)(ak2ak2! 



(a]ci"k2)(a[2aki) 



(99) 



In Eq. (|^ . the first term corresponds to the prod- 
uct of the two single-inclusive distributions, the second 
one gives rise to the Bose-Einstein identical particle cor- 
relation, reflecting the last position, just after particle 
emission. The third term, absent in the tt^tt* case, is 
the responsible for the particle-antiparticle correlation. 
Under analogous conditions as HBT, this last term usu- 
ally vanishes but is non zero if the Hamiltonian of the 
system is of the type H = Hq + Hi, where Hq is the 
free part in the vacuum (final particles) and Hi repre- 
sents the interaction of quasi-particles, resulting in an 
effective shift of their masses. Alternatively, as in the pi- 
oneer work in Ref . (8Qj , this could be similar to having a 
chaotic superposition of coherent states and the density 
matrix containing squeezed states. 

Later, Andreev and Weinerjs^ continued the discus- 
sion. They considered that, in high energy collisions, a 
pion blob (strongly interacting pion system, seen as a 
liquid) is formed and later undergoes a sudden freeze- 
out, where the pionic system (having background and 
excited states) is rapidly converted into free pions. They 
postulate that, at the moment of transition, they can re- 
late in-medium creation and annihilation operators (6^ , b) 
to the corresponding free ones (a^,a), by means of a 
squeezing transformation, with a real squeezing parame- 
ter [r = i \n{Efr/Ein)]- In the same year, Asakawa and 
Csorgo proposed a similar structure to this previous ap- 
proach, but relating in-medium operators to free ones by 
means of a Bogoliubov transformation. They suggested 
that, while the in-medium Hamiltonian, Hi, is diagonal- 
ized in terms of b\ b, the free one is diagonalized by a\d. 
They also proposed to observe hadron mass modification 
in hot medium by means of Back-to-Back Correlations. 
However, the theory was not yet completely correct. 

There were a few more tentative works by the two 
groups but the correct approach was finally written in 
1999, by Asakawa, Csorgo, and Gyulassy^]. 

The formalism developed for the case of an infi- 
nite medium can be summarized as follows. The 
in-medium Hamiltonian is written as H — Hq — 
i/dx dy0(x)(5M2(x - y)0(y), where Hq = 

i / dx ^(/)^ -t- |V(/)p + m^0^j , is the asymptotic (i.e., 

free) Hamiltonian, in the rest frame of matter. The scalar 
field </>(x) in the Hamiltonian H corresponds to quasi 
- particles that propagate with a momentum-dependent 
medium-modified effective mass, which is related to the 
vacuum mass, m, via 

TO^dkl) ^m2-<5M2(|k|). 

The mass-shift is assumed to be limited to long wave- 
length collective modes: (5M^(|k|) ^ if |k| > A^. 



Given such a mass shift, the dispersion relation is mod- 
ified to n^. = uj^ — SAP{\k\), where ilk is the frequency 
of the in-mcdium mode with momentum k. 

The in-medium, thermalized annihilation (creation) 
operator is denoted by 6k (&k), whereas the correspond- 
ing asymptotic operator for the observed quantum with 



four-momentum fc'' = ([Jk,k),a;k 



- k^ (wk > 0) is 



denoted by Ck (Ok). These operators are related by the 

Bogoliubov transformation, i.e., Oki = Cki&ki +slkj6Lk j 
which is equivalent to a squeezing operation. For this 
reason, rk is called mode-dependent squeezing param- 
eter. Note that the relative and the average pair mo- 
mentum coordinates are 5^*2 ~ ~ '^2, qi,2 = ki — k2. 



), andKi,2 = (l/2)(ki -f k2). For short- 



ening the notation, we wrote the squeezed functions as 
Ci,j = cosh[r(i, j , x)] and Sij = sinh[r(?, j, a;)], where 



r{i,j,x) 



-log [{K^^^u,{x))/iK*-<{xW{x))] 



■log 



nk.ix) + Qk, (x) 



(100) 



Also, ni_j is the density distribution, which is taken as 
the Boltzmann limit of the Bose-Einstein distribution, 
i.e., n[*j{x) « cxp{-[/l-*' Uf,{x) ~ ^j.{x)]/T{x)}, where 
the symbol * implies the use of in-medium mass, whereas 
it should be dropped where there is no mass-shift. 

In cases where the boson is its own anti-particle, as 
for tt^tt" or correlations, the full correlation function 
consists of a HBT part (related to the chaotic amplitude, 
Gc(l,2)) together with a BBC portion (related to the 
squeezed amplitude, Gs(l, 2)), as shown below 



C2(ki,k2: 



1 



^ jV2(kl,k2) 
7Vi(ki)7Vi(k2) 

Ge(l,2)|2 ^ |G,(1,2)|2 



G,(1,1)G,(2,2) Gc(l,l)G,(2,2) 



(101) 



being the invariant single-particle and two-particle mo- 
mentum distributions given by 

Gc(l,2) = y/uj^^uj^{al^au^), 
Gs(l,2) = V^kiWk2(akiak2), 
Gc{i,i) = Gc{h,k,)^Ni{k^). (102) 

So far, the derivation considered an infinite system. 
The effects of finite size on BBC will be considered af- 
terwards in the text. In Ref.|8^, they discussed the in- 
fluence of finite emission times, observing that the BBC 
in this case was suppressed when compared to the in- 
stant emission. Including the finite emission time as 
9{t — to)Texp[—r{t~to)], G2(ki,k2) — 1 acquires a multi- 
plicative factor [l + (ajk — w-k)^/r^]- They illustrated the 
maximum BBC, corresponding to G2(k, — k) for the 
meson, obtaining a considerable magnitude for the BBC 
correlation function, in spite of considering the time sup- 
pression. We will sec this more explicitly later. 
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As we discussed in the beginning of this section, the 
BBC is a different type of correlation, discovered for 
boson-antiboson pairs. More recently, Tamas Csorgo, 
Yogiro Hama, Gastao Krein, Prafulla K. Panda, and 
myself demonstrated that a similar correlation existed 
between fermion-antifermion pairs, if the mass of the 
fermions were modified in a thermalized medium. We 
already saw in the Introduction that, regarding the HBT 
effect, identical bosons have an opposite behavior as com- 
pared to identical fermions, as illustrated in Fig. 4, where 
we see that quantum statistics suppresses the probability 
of observing pairs of identical fermions with nearby mo- 
menta, while it enhances such a probability in the case of 
bosons. With respect to BBC, however, we found out a 
very different situation: fermionic BBC are positive and 
similar in strength to bosonic BBC! And, contrary to the 
HBT correlations, the BBC are unlimited. 

In the fermion BBC case, there are expressions similar 
to Eq.jni and (E^, 

iVi(k,) = Wk,(ak,ak,) ; ^i(k,) =Wk.(4.«k,) , (103) 

A^2(ki,k2) oJk^uJi^^ial^al^ai^^ai^^) . (104) 

In the above expressions, < O > denotes the expecta- 
tion value of the operator O in the thermalized medium 
and a\a,d\a are, respectively, creation and annihilation 
operators of the free baryons and antibaryons of mass 



M and lu^ — y + \k^\, which are defined through 
the expansion of the baryon field operator as ^'(x) = 

V EA,A',fe("A,fe«A,fe + «A',-fcaL,_fe)e'''-"; V is the volume 
of the system, ^ and v-^, _^ are the Dirac spinors, 
where the spin projections are A, A' = 1/2, —1/2. The in- 
medium creation and annihilation operators are denoted 
by h\ b, b\ b. While the a-quanta are observed as asymp- 
totic states, the ^-quanta are the ones thermalized in the 
medium. They are related by a fermionic Bogoliubov- 
Valatin transformation. 



aA,k 

4'.-k 





A.k 



here ci = cosri, si = sinri, and 



tan(2ri) 



|ki|AM(ki 



w(ki)2 - MAM(ki) 



(105) 



(106) 



is the fermionic squeezing parameter. Note that in the 
fermionic case, the squeezing parameter is the coefficient 
of sine and cosine functions, differently than the bosonic 
cases in which appeared their hyperbolic counterparts. 
In Eq. (|105|l A is a 2 x 2 matrix with elements A^-^x^ — 
x\^(J ■'kix-^^, where ki = ki/|ki|, x is a Pauli spinor and 
X = —ia'^x- Since r is real in the present case, we drop 
the complex-conjugate notation in what follows. 

In order to evaluate the thermal averages above, the 
system is modelled as a globally thermalized gas o 



quasi-particles (quasi-baryons). In this description, the 
medium effects are taken into account through a self- 
energy function, which, for a spin-^ particle (we will fo- 
cus on proton and anti-proton pairs), under the influence 
of mean fields in a many-body system, can be written as 
E = -|-7°I]° -|-7*E\ In this expression, is a weakly 
momentum-dependent function which, for locally ther- 
malized systems that we are considering, has the role of 
shifting the chemical potential, i.e., yit* = /i — E". The 
vector part is very small and is neglected. The scalar 
part can be written as E** = AM(k). With these ap- 
proximations we describe the system with a momentum- 
dependent in-medium mass, M*(k) = M — AA/(|k|). 

We are mainly interested here in the study of the 
squeezed correlation function, which corresponds to con- 
sidering only the joint contribution of the first and third 
terms on the rhs of Eq. HlOlfl . In the fermionic case 
and for an infinite, homogeneous thermalized medium, 
the correlation BBC part of the correlation function is 
written as 



c: 



(+-) 



{ 



(ki,-ki) = 1+ [l + (2At^k)']"' 

(1 - Wk - nk)^(ckSk)^ 

[c^^k + sl{l - hk)] [clhk + 4(1 - '^k)] 



}, (107) 



where rik - cKp[in^-\,)/T]+i ' ""k - cxp[(Ok+^^^.)/T]+l 
in terms of which the net baryonic density is written as 
Pb = {9/V)J2k ("k - «-k)- In Eq. H107|l we have 
included a more gradual freeze-out by means of a finite 
emission interval, similarly to what was done in Ref [s^jl . 
which has the effect of suppressing the BBC signal. 

For our numerical study of the fermionic back-to- 
back correlations (fBBC), we considered, for simplicity, 
momentum independent in- medium masses, i.e., = 
M — AM. There is no difficulty in considering momen- 
tum dependent self-energies, however, this requires the 
commitment to a specific model and we preferred to leave 
it for a future investigation. In Fig. 25 we show fBBC 
for pp pairs as a function of the in- medium mass M* , for 
three values of the net baryonic density pB- for the nor- 
mal nuclear matter, one tenth of this value and for the 
baryon free region, i.e., pB — 0. We show in the same 
plot results for the bosonic case, bBBC, corresponding 
to (p meson pair, whose mass is close to the proton mass 
and was the example used in Ref.^8^. 

We see from Fig. 25 that fBBC and bBBC are, indeed, 
both positive correlations, with similar shape, and of the 
same order of magnitude. We also observe that fBBC 
is strongly enhanced for decreasing net baryonic density, 
being maximal for « 0, i.e., for approximately equal 
baryon and anti-baryon densities. Fig. 25 also shows that 
fBBC increases with increasing momentum. Besides, we 
see that the shape of fBBC is very sensitive to the shape 
of the freeze-out distribution and, in the limit of a very 
long freeze-out, both fBBC and bBBC would vanish. 
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FIG. 25: The plot on the left shows back-to-back correlations 
of pp (fBBC) and of (/>-meson pairs, for |fc| = 800 MeV/c. 
The plot on the right shows the dependence of the fBBC on 
the in- medium modified proton mass, Af, , and on the net 
baryon density, ps and two typical values of |fc|. In both cases 
T = 140 MeV and At = 2 fm/c. The plots were extracted 
from Ref.lil. 



ing R the radius of the system, T = T{x) the tem- 
perature of the system in each space-time point x, and 
^0 a constant. We assume that the system expands 
with four-dimensional flow velocity u'^ = 7(1, v), where 
V =<u> In the non-relativistic limit, we can write 
7 = (1 + v)^^/^ « 1 + iv^, thus taking into account 
all terms up to O(mv^). We estimate the geometrical 
and dynamical effects on the BBC in the bosonic case, 
considering the in-medium changes for the </)-meson. 




1000 BBC 7.51 




So far, we considered an infinite and homogeneous 
medium, although we know that the systems produced in 
high energy collisions, including the ones at RHIC, have 
finite sizes. Thus, we need to know if the BBC signal 
would survive when more realistic spacial and dynamical 
considerations were made. For pursuing this purpose, we 
studied the effects on the squeezing parameter and on the 
back-to-back correlation of a finite size medium moving 
with collective velocity. For this, we considered a hydro- 
dynamical ensemble, in which the amplitudes Gc and Gs 
in Eq. (|101|l and H102|) are extended to the special form 
derived by Makhlin and Sinyukov [s^ . 



Gc(l,2) 



1,2 



(2^)3 J " 

-I- |s_i,_2p(n_i,_2 + l)}, (108) 
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FIG. 26: The Back-to-Back Correlation is shown as a function 
of the pair average momentum, for ki = — k2 = k, when 
considering the values of the shifted mass around the vacuum 
mass, = 1020 MeV). The squeezed mass was considered 
to change in the entire Gaussian volume, with rms R = 6 fm. 
In the top plots no flow was included, i.e., v =<u> = 0, 
but radial flow with <it>= 0.5 was included in the bottom 
cases. Plots on the left side correspond to the case of particles 
emitted instantly, and those on the right, to the case where 
particles are emitted during a finite time interval, At. 



+ Ci,_2sl2,l(ni,-2 + l)}. (109) 

In Eq.lUnHl and HW^ d^a^'{x) = (fiT.''{x]Tf) F{Tf)dTf 
is the product of the normal-oriented volume element 
depending parametrically on (the freeze-out hyper- 
surface parameter) and on the invariant distribution of 
that parameter F{Tf). We consider two possibilities: i) 
an instant freeze-out, corresponding to F{t) = (5(t — tq); 
ii) an extended freeze-out, with a finite emission inter- 
val, with F{t) = [e{T - To)/At]e-^^~^'>^/^K These cases 
lead, after performing the integration in dr in Eq. H108() 
and H 109(1 with weight {Eij e~*^^''J'^), respectively to: i) 
(lu, + u,) e-*(-'+-.)-o; ii) (w, -t- tOj)[l + [{to, + u;j)At]-\ 
According to the hydrodynamical solution, we can ex- 
press the chemical potential as — — ^j-, be- 



For small mass shifts, i.e. ^ 1, the flow ef- 

fects on the squeezing parameter are of fourth order, i.e., 
^Kj^n_^^ergy'j (^s^y ^ consequcncc, the flow ef- 
fects on Ti j can be neglected, and the factor Cij and Sij 
become flow independent, although they could still de- 
pend on the coordinate r through the shifted mass, m* 
(e.g., through T(x), as in hydrodynamics), which is not 
considered here. 

For the sake of simplicity, trying to keep the results as 
much analytical as possible (for details, see |1^), we make 
the hypothesis that the mass-shift is independent on the 
position within the fireball. We further assume that this 
last one has a sharp surface, i.e., Sm = on the surface, 
and also the density vanishes outside the system volume. 
The spatial integration in Eq. p()8|l and H109|l extends 
over the region where the mass-shift is non-vanishing, 
which is not inflnitcly large. For instance, I in relativistic 
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heavy ion collisions is a finite region y w i?"^ « (5 — 10)'^ 
fm'^. We should keep in mind that the vacuum term 
in the integrand vanishes outside the mass-shift region, 
since it is proportional to s^j , which vanishes in that 
region. On the other hand, the terms proportional to 
n'fj (x) are finite. Being so, we can extend the integration 
in Eq. H108|) and (|109|l to infinity and, without much loss 
of generality, we can choose for V a Gaussian profile, 
exp[-rV(2i?2)]. 

We studied the flow effects in two cases. The first one 
considers the mass-shift as occurring in the entire volume 
of the system, for which the parameter R, in the exponen- 
tial, represents the cross-sectional area of the Gaussian 
profile with rms R . The other case corresponds to con- 
sidering the shift of mass in a smaller volume, associated 
to a certain Rg < R. We show in Fig. EHlthe results for 
the first case only, where we adopted i? = 6 fm. The BBC 
signals for the mass shift in the partial volume, although 
not shown, are very similar to the ones in Fig. 26, except 
for the maximum strength of the bBBC signal, which is 
smaller than seen in the plots of Fig. 26, emphasizing 
that the effect is directly proportional to the size of the 
region where the mass-shift occurs. By comparing the 
top panels (no flow) with the bottom ones (<u>= 0.5) in 
Fig. 26, we see that the flow has a suppression effect on 
the bBBC strength in the high momentum region. How- 
ever, the presence of flow surprisingly enhances the BBC 
signal in the low momentum region (i.e., for |k| ^ 1000 
MeV) , as compared to the no-flow case on top panel, sug- 
gesting this as the region to focus the search for the BBC 
effect. In any case, we conclude with enthusiasm that the 
BBC signal survives more realistic conditions, like finite 
sizes, time spread, and flow. This strongly encourages 
us to analyze even more realistic systems than the ones 
discussed here, and to optimize the way we should look 
for this interesting, although not yet observed, effect. 



BRIEF CONCLUDING REMARKS 

I would like to finalize by emphasizing that, although 
a lot was written above about the evolution and contri- 
butions related to HBT, since its discovery about half a 
century ago, the present text merely covers a tiny fraction 
of what has been produced on this subject. An increas- 
ing number of people over these five decades have given 
countelss contributions to the field. Unfortunately, the 
lack of space would not allow me to mention and discuss 
them all. Finally, I would like to dedicate this review to 
the memory of my father, who passed away this year, as 
well as to Robert Hanbury-Brown and Richard Q. Twiss, 
on celebrating the 50*'* anniversary of their first publica- 
tion on this fabulous method and discovery. 

I would like to express my gratitude to Yogiro Hama 
and Takeshi Kodama for devising, 15 years ago, and orga- 
nizing ever since, this active group working on hadronic 



and nuclear interactions at high energies. 
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